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I.  INTRODUCTION 

The  objectives  of  this  contract  are:  (1)  to  relate  the  6300  A emissions 
as  observed  by  the  ISIS-2  scanning  photometer  to  the  causative  electron  fluxes 
which  are  simultaneously  observed  by  the  ISIS-2  Soft  Particle  Spectrometer. 

The  energy  sub-region  of  the  electron  flux  which  is  primarily  responsible  for 
the  6300  A emissions  will  be  delineated  to  the  accuracy  available  with  the 
ISIS-2  instrumentation.  It  is  hoped  that  an  empirical  relation  can  be  derived 
that  relates  electron  energy  flux  to  6300  A flux.  (2)  To  compare  the  particle 
fluxes  and  6300  A emissions  (and  the  conclusions  derived  in  Item  1)  with  simul' 
taneous  electron  density  profiles  obtained  by  the  ISIS-2  topside  sounder.  The 
basic  goal  of  such  a comparison  will  be  to  ascertain  if  it  is  possible  to 
infer  gross  characteristics  of  the  ionosphere  (e.g.  fQF2*  P^^sma  trough  boun- 
daries) from  6300  A and  particle  flux  measurements.  (3)  To  relate,  where 
possible,  the  results  obtained  in  items  1 and  2 to  the  dynamical  effects  pro- 
duced by  substorms.  Specifically  we  will  evaluate  if  it  is  possible  to  infer 
the  motion  of  boundaries  such  as  the  poleward  edge  of  the  low-altitude  plasma 
trough  from  indirect  measurements.  (4)  To  analyze  the  simultaneous  data  ob- 
tained with  ISIS-2,  the  akol  Airborne  Ionospheric  Laboratory,  and  the  Defense 
Meteorological  Satellite  System.  Again  the  goal  will  be  to  ascertain  if  iono- 
spheric parameters  can  be  determined  by  indirect  means. 

This  report  represents  the  work  accomplished  during  the  second  twelve 
months  of  a 39  month  contract.  Most  of  the  work  accomplished  during  the 
first  24  months  was  preparatory  in  nature.  However,  items  1 & 3 of  the 
objectives  have  been  brought  to  a reasonable  conclusion. 
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II.  ACCOMPLISHMENTS 

Item  1 of  the  contract  states  that  one  objective  is  "to  relate  the 
6300  A emissions  as  observed  by  the  ISIS-2  scanning  photometer  to  the 
causative  electron  fluxes  which  are  simultaneously  observed  by  the  ISIS-2 
Soft  Particle  Spectrometer.  The  energy  sub-region  of  the  electron  flux 
which  is  primarily  responsible  for  the  6300  A emissions  will  be  delineated 
to  the  accuracy  available  with  the  ISIS-2  instrumentation.  It  is  hoped 
that  an  empirical  relation  can  be  derived  that  relates  electron  energy 
flux  to  6300  A flux." 

Originally  we  had  envisioned  the  necessity  of  determining,  indepen- 
dently, the  relationship  between  several  energy  bands  and  the  6300  A 
fluxes.  This  would  have  required  defining  regions  where  one  energy  band 
would  dominate  in  the  production  of  6300  A.  However,  this  turned  out  to 
be  unneccessary,  as  will  be  described  below. 

In  order  to  expedite  the  analysis  of  a large  number  of  passes,  UTD 
and  York  University  developed  programs  to  merge  the  data  from  the  SPS, 

ASP  (5577  A and  3914  A),  and  the  RLP  (6300  A)  instruments.  The  energy 
flux  within  the  local  loss  cone  was  evaluated  from  the  SPS  data  and 
computer  plotted  along  with  the  three  optical  lines.  Figure  1 shows  an 
example  of  this  plot.  Four  energy  bands  are  selectable  at  one  time. 

Ten  passes  were  chosen  for  an  initial  evaluation  of  the  technique. 
Inspection  of  these  passes  revealed  the  comparison  problem  to  be  much 
easier  than  originally  envisioned.  It  was  discovered  that  to  first 
order  > 300  electron  fluxes  produce  a negligible  amount  of  6300  A 
emission.  This  can  be  easily  seen  in  Figure  1.  Between  63**  and  68° 
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invariant  latitude  there  is  a two  orders  of  magnitude  increase  and  decrease 
in  the  c 'eater  than  300  eV  energy  flux  with  no  corresponding  change  in  the 
6300  A emission  rate;  5577  and  3914  A do  track  the  > 300  eV  flux  increase. 
All  ten  passes  examined  exhibit  this  relationship.  Thus  only  the  < 300  eV 
fluxes  need  be  considered  in  any  relationship  between  electrons  and  6300  A 
emission. 

At  present  we  are  expanding  the  number  of  passes  to  approximately  50. 
An  empirical  relationship  between  < 300  eV  electron  energy  fluxes  and 
6300  A emission  rates  will  be  developed  from  this  data  set  during  the 
next  year,  thus  satisfying  item  one  of  the  contract. 

As  discussed  in  the  previous  yearly  report  a relationship  between  the 
IMF,  and  AE  and  the  equatorward  cessation  of  >100  eV  electron  fluxes  has 
been  derived.  (The  paper  describing  these  results  will  appear  in 
shortly.)  Combining  the  results  from  this  paper  and  the  work  described 
above  the  equatorward  boundary  of  auroral  6300  A can  be  related  to  the 
external  IMF  and  AE,  thus  satisfying  Item  3 of  this  contract. 

During  the  past  year  contour  plots  of  6300  A,  5577  A,  and  3914  A 
emission  rates  have  been  submitted  to  the  World  Data  Center.  Ninety-five 
ISIS-2  passes  were  included  in  this  submission. 
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Abstract . Results  are  presented  from  low-energy 
plasma  analyzers  (12  eV  to  12  keV)  carried  on 
two  rockets  launched  into  the  dayside  cleft 
during  January  1975.  We  conclude  that  (1) 
atmospheric  interaction  becomes  Important  for 
<l-keV  electrons  at  approximately  250  km,  (2) 
characteristics  of  particles  In  'Inverted  V's' 
observed  in  the  afternoon  cleft  are  consistent 
with  their  Interpretation  as  being  due  to 
parallel  electric  field  acceleration  from  a 
constant  source  population,  and  (3) 
magnetospheric  'energetic'  (>2  keV)  electrons 
Intermingle  with  'magnetosheathlike'  plasma  in 
the  cleft. 

Introduction 

The  past  two  years  have  marked  a dramatic 
upsurge  in  sounding  rocket  launches  Into  the 
low-altitude  cleft  region.  Prior  to  1974  only 
four  instrumented  rocket  payloads  {Maynard  and 
Johns  tone . 1974;  Ledley  and  Farthing.  1974] 
and  one  thermite  barium  release  [Mlkkelsen  and 
Jorgensen,  1974]  had  been  launched  into  or  In 
the  vicinity  of  the  cleft.  The  Instrumented 
payloads  were  launched  from  Hall  Beach, 
Northwest  Territories  (A*  79.7*),  on  March  15, 
18,  19,  and  22,  1971,  and  Che  barium  release 
from  Sondre  Stromfjord,  Greenland  (A-  74.5°),  on 
December  10,  1972. 

In  contrast,  from  June  1974  to  November  1975, 
15  rockets  were  fired  in  the  cleft  region  from 
launch  sites  at  Cape  Parry,  Northwest 
Territories,  Canada,  and  Sondre  Stromfjord, 
Greenland.  Ungstrup  et  al.  [1975a]  reported 
results  from  two  rockets  launched  from  Greenland 
in  July  1974;  these  payloads  were  more 
comprehensive  than  their  predecessors  as  they 
included  particle  spectrometers,  dc  electric  and 
magnetic  probes,  and  thermal  plasma  sensors. 
During  December  1974  through  January  1975,  two 
major  c.-impalgns,  totaling  11  rockets,  were 
executed  in  Canada  (Cape  Parry,  Northwst 
Territories, A • 74.5°)  and  Greenland  (Sondre 
Stromfjord).  These  launchings  li.,’"ded 
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extensive  plasma  and  field  instrumentation, 
optical  sensors,  shaped  charge  barlim\ 
injections,  and  lithium  releases  [Haerendel . 
1975;  Torbet  et  al..  1975;  Ungstrup  et  al., 
1975b;  Temerln  et  al..  1975;  Shepherd  et  al.. 
1975,  1976;  Wlnnlngham  et  al.,  1975a;  Wescott  et 
al. , 1975;  Stenbaek-Nlelsen  et  al..  1975; 
Jeffries  et  al.,  1975;  Torbet  and  Carlson, 
1975]. 

Another  series  consisting  of  two  rocket 
launches  was  completed  during  November  1975  at 
Cape  Parry.  These  vehicles  Included  particle 
detectors,  barium  lined  shaped  charges,  plasma 
drift  detectors,  and  energetic  particle  sensors. 

In  this  paper  we  «d.ll  present  results  from 
the  soft  particle  spectrometer  (SPS)  [Helkklla 
et  al..  1970]  instrisnents  carried  as  an 
environmental  monitor  on  Tordo  Uno  and  Dos.  The 
Tordo  Uno  and  Dos  programs  %fere  designed  by  the 
Los  Alamos  Scientific  Laboratory  and  the 
University  of  Alaska  Geophysical  Institute,  and 
the  rockets  were  launched  from  Cape  Parry  by  the 
Sandla  Laboratories  on  January  6 and  11,  1975, 
2349:02  and  0025:02  UT,  respectively. 


Instrumentation 

The  primary  objective  of  the  Tordo  campaign 
was  the  injection  of  barium  ions  Into  the 
low-altitude  cleft  (Initial  results  from  the 
shaped  charge  barium  releases  are  described  In 
the  paper  by  Jeffries  et  al.  [1975]).  In 
addition  to  the  shaped  charge,  a low-  energy  (12 
eV  to  12  keV)  plasma  analyzer  (SPS)  similar  to 
the  one  on  the  Isis  1 and  2 satellites  [Heikklla 
et  al. . 1970]  was  carried  as  a monitor  to  locate 
the  detonation  point  relative  to  the  cleft.  For 
proper  barium  Injection  payloads  were  actively 
stabilized  In  relation  to  the  local  magnetic 
field,  the  result  being  a unique  SPS  pitch  angle 
of  45  . The  SPS  concurrently  measured  the 
differential  energy  spectrum  of  both  electrons 
and  positive  Ions  In  15  logarithmically  spaced 
steps  from  12  eV  to  12  keV  every  3.2  s.  In 
addition,  there  was  a sixteenth  step  at  zero 
energy  to  monitor  the  background  count  rate. 
Each  energy  level  was  held  for  0.2  s,  and  the 
count  rate  was  sampled  eight  times  during  this 
period.  Multiple  sampling  was  employed  in  order 
to  evaluate  possible  temporal/spatial  aliasing 
occurring  during  the  3.2  s required  to  complete 
a stepping  sequence. 

The  energy  band  pass  of  the  SPS  was  set  at 
32T  (being  a divergent  electrostatic  lens,  the 
SPS  resolution  is  determined  by  a field  stop), 
and  the  geometric  factor  (not  Including  the 
energy  band  pass)  was  8.6  x lo”^  cm  sr.  Ir. 
order  to  match  the  energy  band  passes  at  full 
width  at  half  maximum  (FWHM),  22  logarthmlcal Iv 
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spaced  steps  would  be  required  instead  of  the  13 
used  on  these  flights. 

The  energy-analyzed  particles  were  detected 
by  lA-stage  RCA  discrete  dynode  electron 
multipliers  [Bunting,  et  al.,  19721,  and  the 
subsequent  pulses  were  .'upllfled  by  high  gain 
bandwidth  amplifiers  [Smith,  1972] . 
Postacceleration  was  employed  to  Increase  the 
detection  efficiency  of  the  multipliers  at  low 
energies.  The  pulse  train  from  each  amplifier 
was  log2  coiapressed  and  stored  for  subsequent 
telemetry  to  the  ground. 

Magnetic  Conditions 

Tordo  Uno  was  launched  on  January  6,  1975,  at 
2349:02  UT  (approximately  1330  magnetic  local 
time).  The  Kp  value  for  the  3-hour  period 
including  Che  flight  was  6+  and  for  the 
following  period  was  8-.  Figure  la  gives  All  and 
AI.  for  the  6-hour  period  bracketing  the  Tordo 
Uno  flight.  It  should  be  noted  that  the  AE 
Index  shown  Is  based  on  the  standard  11  auroral 
zone  stations  plus  data  from  an  additional  8 
stations  in  Greenland  (courtesy  of  T.  S. 
Jorgensen)  and  Scandanavla  (courtesy  of  R.  S. 
Pelllncn)  . Inspection  of  Figure,  la  Indicates 
that  the  magnetosphere  was  In  an  extreme  state 
of  agitation  during  this  6-hour  period.  During 


Fig.  la.  AU  and  AL  Indices  for  January  6, 
1973.  The  auroral  Index  Is  constructed  of 
the  standard  11  auroral  zone  stations  plus 
an  additional  8 stations  on  the  nlghtslde 
in  Greenland  and  Scandanavla. 


the  flight,  AE  was  as  large  as  750Y.  Solar  wind 
magnetic  field  data  were  not  available  for  this 
period. 

Tordo  Dos  was  launched  on  January  11,  1973, 
at  0023:02  UT  (approximately  1400  magnetic  local 
time).  The  two  days  encompassing  the  Tordo  Dos 
flight  were  relatively  quiet;  January  10  was  a Q 
day  and  January  11  was  a QQ  day.  The  Up  value 
for  the  last  3-hours  of  January  10  was  1+,  and 
for  the  first  3-hours  of  January  11  the  Kp  value 
was  2.  Examination  of  Figure  lb  shows  that 
magnetic  activity  was  weaker  and  more 
representative  of  classical  substorm  morphology 
during  the  6-hour  period  bracketing  Tordo  Dos 
than  during  Tordo  Uno,  The  flight  of  Tordo  Dos 
occurred  during  the  recovery  phase  of  a small 
high-latltude  substorm  (maximum  AE  ■ 250Y). 

Interplanetary  magnetic  field  data  from  Imp  8 
were  available  during  the  Tordo  Dos  flight  and 
are  presented  In  Figure  2.  The  IMF  went 
southward  at  2216  UT  on  January  10  and  remained 
southward  or  close  to  zero  until  0102  on  January 
11,  when  it  became  positive  again.  A more 
detailed  description  of  the  motion  of  the  cleft 
and  Its  relation  to  the  IMF  during  January  10 
and  11  can  be  found  In  the  work  of  Stiles  et  al. 
[19771  . 

Experimental  Results 

Tordo  Uno.  As  was  mentioned  In  the  previous 
section,  Tordo  Uno  was  launched  during  a period 
of  prolonged  magnetic  activity.  Examination  of 
the  bottomslde  sounder  records  from  Cape  Parry 
and  Sachs  Harbour  (see  Jeffries  et  al.  [1975] 
and  Stiles  et  al.  [1977]  for  a more  complete 
description)  Indicated  that  the  equatorward  edge 
of  the  cleft  was  south  of  Cape  Parry.  The 
general  characteristics  of  the  lonograms 
Indicated  precipitation  of  particles  at  both 
Parry  and  Sachs.  As  a result  of  the  above 
considerations  a decision  was  made  to  launch 
westward  In  order  to  avoid  possible  exiting  the 
precipitation  region  northward  Into  the  polar 
cap.  Figure  3 shows  the  projection  of  the 
rocket  trajectory  to  100  km  along  the  earth's 
magnetic  field.  From  Figure  3 It  can  be  seen 
that  the  rocket  traveled  from  approximately  75° 
to  76.4°  Invariant  latitude. 

A synoptic  view  of  the  electron  data  from 
Tordo  Uno  is  presented  In  the  energy-time 
spectrogram  format  In  Figure  4.  This  format  Is 
essentially  the  same  as  that  used  for  Isis  1 and 
2 spectrograms  [see  Winningham  et  al.,  1975b], 
with  the  addition  of  the  altitude  profile  (solid 
line)  to  the  energy  flux  panel  and  the  launch 
universal  date  and  time  at  top  left.  The  Z axis 
(grey  scale)  of  the  spectrogram  (top  panel)  has 
been  scaled  to  give  the  count  per  11. T1  ms 
Instead  of  the  actual  acumulatlon  period.  This 
was  done  In  order  to  have  the  same  relationship, 
between  count/accumulation  period  and  grey 
scale,  as  that  In  the  Isis  1 and  2 spectrograms. 

One  of  the  most  striking  features  observed  in 
the  electron  spectrogram  Is  the  absence  of  major 
spatlal/tenporal  structure.  This  was  not  the 
case  for  Tordo  Dos,  as  will  be  seen  In  the 
following  section.  Throughout  the  Tordo  Uno 
flight  the  count/accumulation  period 
(differential  energy  flux)  maximized  at 
approximately  100-200  eV. 
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Fig.  lb.  Same  as  Figure  la  but  for  January  11,  1975. 


A 6300-A  scanning  photometer  belonging  to  G. 
C.  Shepherd,  York  University,  and  an  all-sky 
camera  were  operated  at  Sachs  Harbour  during  the 
Tordo  campaign.  The  photometer  indicated  a 
uniform  glow  of  about  300  R,  increasing  to  1 k R 
toward  the  south  of  Cape  Parry  during  the  flight 
of  Tordo  Uno.  Also  the  all-sky  camera  data 


Independent  (Figure  7).  At  2353:26.4  UT  an 
abrupt  increase  in  magnitude  was  observed  in  all 
energy  channels  below  1021  eV.  The  only 
spectral  change  occurring  was  a slight 
enhancement  in  the  lowest  energy  channel  (12 
eV) . At  2354:30.4  UT  an  abrupt  decrease  was 
observed,  after  which  time  the  flux  remained  at 


showed  a lack  of  any  discrete  aurora, 

A more  quantitative  presentation  of  the  data 
is  given  in  Figures  5,  6,  and  7.  Figure  5 
presents  electron  differential  flux  versus 
altitude,  Figure  6 a representative  electron 
spectrum  near  apogee,  and  Figure  7 differential 
flux  versus  time.  Inspection  of  Figure  5 and  7 
reveals,  as  did  inspection  of  Figure  4,  the 
absence  of  dramatic  acceleration  events  (see 
Figure  8 for  comparison) . 

Throughout  the  flight  the  electron  spectrum 
can  be  divided  into  two  relatively  distinct 
regions,  above  and  below  approximately  2 keV. 
Below  2 keV  the  spectrum  appears  to  be  comprised 
of  a low-energy  'primary'  plus  'secondary' 
population,  the  primary  portion  of  the  spectrum 
being  broader  than  a Maxwellian  curve.  For 
reference,  a Maxwellian  curve  has  been  included 
in  Figure  6 and  has  been  normalized  to  the  87-eV 
experimental  point  (most  of  the  spectra  around 
apogee  exhibit  a relative  maximum  at  87  eV) . 
Also  included  in  Figure  6 are  representative 
near-concurrent  plasma  man'le  and  raagnetosheath 
spectra  measured  by  Hawkeye  (courtesy  of  L,  A, 
Frank).  The  Hawkeye  data  will  be  considered 
again  in  the  discussion  section. 

The  electron  spectrum  above  2 keV  was 
statistically  poor  but  well  above  the  noise 
level.  Typical  counting  rates  were  less  than 
100/s.  The  spectrum  presented  in  Figure  6 is  an 
average  of  26  spectra  (83.2  s)  beginning  at 
2350:36.8  and  ending  at  2352:00  UT.  The 
deviation  of  each  individual  spectrum  from  the 
average  was  minor.  This  high-energy  portion  of 
the  electron  spectra  Increased  and  decreased 
throughout  the  flight  (see  Figure  7a) 
Independently  of  the  lower  energy  population. 

Examination  of  Figures  Sa  and  5b  reveals  a 
sharp  increase  of  flux  with  altitude  between  180 
and  250  km.  This  altitude  related  variation  was 
energy  dependent,  loirer  energies  displaying  the 
largest  Increase  with  increasing  altitude. 
Above  250  km  (2350:56  UT)  the  dependence  was 
mainly  spatial  and/or  temporal,  not  altitudinal. 
From  2350:56  to  2353:26.4  UT  there  was  a steady 
increase  in  the  electron  flux  at  energies 
between  87  and  1021  eV  that  was  energy 


values  close  to  the  preincrease  level  except  for 
three  narrow  bursts. 

At  energies  below  87  eV  an  initial  flux 
increase  with  altitude  was  also  observed. 
However,  at  later  times  the  flux  profiles  were 
not  Identical  to  those  at  higher  energies 
discussed  earlier.  The  initial  slow  increase 
was  observed  in  the  52.8-eV  channel  but  not  in 
the  30. 8-,  18. 9-,  and  11.7-eV  channels.  The 
30.8-eV  sample  exhibited  a relatively  flat 
profile  with  altitude  prior  to  the  Increase.  At 
18.9  and  11.7  eV  a relative  maximum  was 
observed,  the  18.9-eV  sample  leading  the  11.7-eV 
sample  in  either  apace  or  time  (2351:30  and 
2352:00  UT,  respectively).  As  in  the  higher 
energy  channels,  an  abrupt  Increase  and  decrease 
in  flux  was  observed  after  2353:26.4  UT. 

Throughout  the  flight  positive  ion  fluxes 
were  weak.  At  apogee  the  peak  ion  flux  occurred 
in  the  few  hundred  electron  volt  range  with  flux 
levels  of  approximately  lO"*  cn“‘  sr”  s”  . 

Tordo  Dos.  lonosonde  returns  obtained  prior 
to  launch  indicated  that  the  cleft  was  poleward 
of  (^pe  Parry  and  close  to,  but  poleward  of, 
Sachs  Harbour.  (The  lonosonde  returns  for  this 
day  are  described  in  detail  by  Stiles  et  al. 
(19771.)  The  decision  had  been  reached  that 
Tordo  Dos  would  be  launched  across  the 
equatorward  boundary  of  the  cleft.  In  order  to 
accomplish  this  goal  the  launch  azimuth  was  set 
as  close  to  magnetic  north  as  range  safety 
requirements  would  allow  (Figure  3). 

From  high  voltage  turnon  at  0026:27.6  UT  (180 
km)  to  0028:13.2  UT  (390  km),  downcoming 

IMP  a 

INTERPLANETARY  MAGNETIC  FIELD 


UTZlOO  2200  2J00  2400  OiOO  0200 


10-11  JANUARY  I9T6 

Fig.  2,  Interplanetary  magnetic  field  Z 
component  (solar  magnetospheric)  from  Imp  8 
for  January  10-11,  1975. 
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Fig,  3.  Tordo  Uno  and  Dos  tra.iectorles  projected  Co  100  km  along  the  earth's 
magnetic  field. 


atmospheric  photoelectrons  and  very  weak 
high-energy  electrons  were  observed  (Figure  8). 
^ maximum  In  the  photoelectron  flux  was  observed 
between  approximately  300  and  340  km.  Mantas 
and  Bowhlll  [1975]  predict  a maximum  In  the 
downcoming  photoelectron  flux  at  approximately 
290  km  for  a 9o"  solar  zenith  angle.  Our  peak 
was  above  this  altitude,  but  the  zenith  angle 
(100  ) was  also  larger.  In  addition,  electron 
density  fluctuations,  causing  the  observed 
spread  F,  would  result  In  fluctuations  In  the 
atmospheric  photoelectron  flux  and  a smearing  of 
the  altitudinal  profile  (G.  P.  Mantas,  private 
communication,  1976). 

F.xamlnatlon  of  the  SPS  data  Indicates  that 
the  cleft  equatorward  boundary  was  crossed  at 
002H;13.2  UT  at  an  Invariant  latitude  of  76.4° 
(Figure  9).  From  0028il3.2  to  0029:10.8  UT  the 
electron  spectrum  peaked  at  approximately 
100-200  eV  and  varied  In  magnitude.  The 
spectral  shape  Is  similar  to  shapes  reported  In 
the  literature  for  the  cleft  IWlnnlngham  and 
Helkkl la.  1974).  Subsequent  to  this  time,  two 
'Inverted  V's'  were  observed,  one  centered  at 
0029:26.8  and  the  other  at  00:31  UT. 
Morphologically,  inverted  V's  are  regions  In 
energy-time  spectrograms  where  the  maximum  count 


rate  (region  of  maximum  darkness  or  lightness, 
depending  of  film  polarity)  moves  first  to 
higher  energy  then  to  lower  energy  with 
advancing  time,  producing  an  Inverted  V shape. 
Between  the  two  Inverted  V's  the  electron 
spectra  were  either  peaked  In  the  100-  to  200-eV 
range  or  were  due  to  atmospheric  photoelectrons. 

Inspection  of  each  electron  spectrum  during 
passage  through  the  Inverted  V's  Indicated  that 
on  a gross  scale  the  peak  In  the  electron 
distribution  function  (particles  cm“^  sr”^  eV'^ 
s ) was  increasing  and  then  decreasing  in  both 
energy  and  magnitude.  In  addition,  several 
relative  maxima  were  observed  (see  Figure  12), 
The  spectrum  In  Figure  10a  was  obtained  Just 
prior  to  the  first  'energized'  spectrum  In  the 
second  inverted  V.  It  Is  well  fitted  by  a 
Maxwellian  distribution  (dashed  curved)  with  a 
'temperature'  of  1.7  x 10°  °K  and  a density  of 
0.5  cm”  (asainlng  Isotropy).  The  solid  circles 
represent  a correction  for  atmospheric 
photoelectrons  utilizing  the  spectrum  measured 
at  0030:11.6  UT.  Spectra  In  Figure  10b  are 
representative  of  the  peaked  distributions  In 
the  second  Inverted  V. 

All-sky  camera  data  were  obtained  frcxn  Sachs 
Harbour,  which  at  the  time  of  Tordo  Dos  had  just 
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Fig.  4.  Electron  energy-time  spectrogram  for  Tordo  Uno. 


entered  astronomical  tvd.llght  (12°  solar 
depression  angle).  In  the  all-sky  camera  data, 
faint  auroras  are  visible  toward  the  north,  and 
their  positions  are  plotted  in  Figure  9 together 
with  those  of  the  two  Inverted  V events  observed 
by  the  rocket  particle  detectors  and  cleft 
returns  deduced  from  the  ionosonde  data.  All 
the  data  shown  have  been  mapped  to  an  altitude 
of  160  km.  The  height  of  160  km  used  for 
calculating  the  position  of  the  auroras  was 
selected  on  the  basis  of  the  energy  spectra 
(Figure  8)  and  luminosity  profiles  derived  by 
Rees  [1963]. 

As  can  be  seen  in  Figure  9,  auroras  were 
observed  only  in  the  vicinity  of  the  second 
Inverted  V event.  No  auroral  form  associated 
with  the  first  event  could  be  discerned  in  the 
data.  This  at  first  appeared  strange 
considering  the  very  similar  energies  and  energy 
fluxes  in  the  two  events.  As  will  be  discussed 
below,  the  explanation  is  to  be  found  in  a 
combination  of  viewing  angles  from  Sachs  Harbour 
and  the  latitudinal  width  of  the  two 
precipitation  events. 

From  Figure  8 we  find  a characteristic  energy 
In  the  two  events  of  1 keV  and  an^energy  flux  of 
about  1 erg/cm  s sr  or  3 ergs/cm s (assuming 
Isotropy).  For  a 1-keV  characteristic  energy, 
Rees  and  Luckey  (1974),  calculate  a column 
emission  rate  of  0.16  kR/erg/cm^  s in  4278  X and 
a corresponding  emission  rate  ratio  5577/4278  X 
of  6.  Thus  the  estimated  brightness  of  the 
form,  if  overhead,  would  be  3 x 0.16  x 6 • 3 kR 


C 

ill  5577A.  The  forms  are  not  overhead  but  are 
off  to  the  northeast  at  an  elevation  angle  of 
about  45°  and  35°,  respectively,  and  the  SPS 
data  (Figure  8)  indicate  the  latitudinal  width 
of  the  precipitation  in  the  second  event  to  be  3 
times  that  of  the  first.  Romick  and  Be  Ion 
[1967]  have  calculated  the  brightness  of  an 
aurora  for  various  values  of  width,  position 
relative  to  the  observer,  and  height  luminosity 
profile.  They  find  that  the  brightness  of  an 
arc  0.3°  wide  at  an  elevation  angle  of  35°  (very 
similar  to  the  values  for  the  second  Inverted  V 
event)  would  be  0,35  of  the  overhead  value. 
Further,  they  find  that  auroras  well  away  from 
the  zenith  have  a brightness  roughly 
proportional  to  their  thickness,  all  other 
things  being  equal.  With  these  corrections  for 
aspect  and  width  applied  the  auroras  associated 
with  the  two  Inverted  V events  would  appear  with 
a brightness  of  0.3  and  1 kR  with  respect  to  an 
observer  at  Sachs  Harbour. 

Under  ideal  conditions  the  35-mm  f 1.2 
all-sky  camera  used  at  Sachs  Harbour  can  resolve 
forms  of  slightly  less  than  1 kR  In  5577  X 
(Romick  and  Brown.  19711.  Thus  the  aurora 
associated  with  the  first  Inverted  V event  (0.3 
kR)  should  not  be  visible  In  the  all-sky  data 
even  under  Ideal  conditions,  and  given  the 
twilight  condition  during  the  experiment,  the 
aurora  (1  kR)  associated  with  the  second 
Inverted  V should  only  be  marginally  detectable, 
which  is  in  good  agreement  with  the  actual  data. 

It  should  be  noted  that  the  6300-X  photometer 
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Fig.  5.  Differential  number  flux  versus  altitude  for  Tordo  Uno. 

observed,  in  addition  to  the  auroral  form  also  cleft.  Comparison  of  Tordo  Uno  positive  ion 

detected  by  the  all-sky  camera,  a faint  arc  just  fluxes  vd.th  Isis  results  indicate  a position  in 

north  of  Sachs  Harbour  at  an  invariant  latitude  the  poleward  portion  of  the  cleft.  This 

of  76.4°.  The  photometer  was  scanning  conclusion  is  based  on  two  characteristics  of 

geographic  north-south  and  thus  was  pointing  the  positive  ions,  namely,  their  low  intensity 

away  from  the  rocket  trajectory,  and  it  is  and  low  average  energy.  Typical  results  from 

uncertain  what  relevance  these  measurements  have  Isis  show  the  hardest  and  most  intense  fluxes  at 

to  the  observations  along  the  rocket  trajectory.  the  equatorward  edge  of  the  cleft  with  a 

A strong  oblique  E layer  return  is  seen  in  subsequent  poleward  softening  and  decrease  in 

the  Cape  Parry  lonograms,  and  several  returns  Intensity. 

are  also  seen  in  the  Sachs  Harbour  lonograms.  Another  possible  Interpretation  of  the 

By  assuming  that  the  returns  come  from  results  would  be  that  the  observed  electrons  are 

ionospheric  structures  aligned  along  constant  polar  rain  electrons  IWinningham  and  Heikklla. 

magnetic  latitude  a triangulation  can  be  made  to  1974).  However,  several  facts  argue  against 

determine  the  location  of  the  irregularity.  By  this  interpretation.  The  observed  density  of 

triangulating  on  the  strongest  oblique  E layer  ^he  primary  electrons  is  as  high  as 
return  its  position  was  found  to  lie  within  the  approximately  5 cn"^.  Winningham  and  Heikklla* s 

region  of  the  second  Inverted  V event  as  shown  [1974)  results  show  the  polar  rain  fluxes  to  be 

in  Figure  9.  only  a few  percent  of  the  cleft  fluxes.  If 

From  Figure  9 it  can  be  seen  that  the  cleft  these  densities  were  measured  in  the  polar  cap, 

as  measured  by  the  SPS  was  Indeed  poleward  of  cleft  densities  would  be  >100  cm"^  and  by 

Sachs  Harbour  by  a few  tenths  of  a degree.  The  Inference,  magnetosheath  densities  would  be  at 

position  of  the  cleft  as  determined  by  the  SPS  least  as  large.  Such  cleft  densities  are  rarely 

(76.4  ) Is  coincident  with  the  aforementioned  observed.  An  even  more  telling  argument  comes 

weak  red  arc  even  though  they  are  separated  by  concurrent  measurements  of  solar  wind 

approximately  4 of  longitude.  densities  by  Imp  7 (courtesy  of  W.  C.  Feldman) 

As  In  Tordo  Uno,  the  proton  fluxes  were  very  apj  magnetosheath  densities  by  Hawkeye  (Figure 

weak.  It  should  be  noted  that  the  SPS  as  used  (courtesy  of  L.  A.  Frank).  Both  were 

in  Tordo  Uno  and  Dos  was  approximately  20  times  approximately  10  cm"^  during  the  period  of  Tordo 

more  sensitive  than  were  the  Isis  SPS's.  (The  Uno. 

main  Increase  in  sensitivity  is  due  to  the 

increased  accumulation  period.)  Comparison  with  the  Hawkeye  results  (Figure 

6)  indicates  a factor  of  approximately  2 
Discussion  and  Conclusions  decrease  in  density  between  the  magnetosheath 

and  the  poleward  portion  of  the  cleft.  The 
The  exact  location  of  the  Tordo  Uno  results  of  Akasofu  et  al.  (1973)  demonstrate 

measurements  relative  to  the  cleft  boundaries  is  that  the  poleward  portion  of  the  cleft  maps  to 

difficult  to  establish.  Only  one  detail  can  be  the  inner  edge  of  the  magnetospheric  boundary 

stated  with  certainty,  that  Is,  Tordo  Uno  was  layer  (see  their  Figure  11).  Additionally 

poleward  of  the  equatorward  boundary  of  the  Akasofu  et  al.  show  the  boundary  layer 
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Fig.  6.  Electron  spectra  from  Tordo  Uno  and  Hawlceye  for  January  6,  1975. 


densities  to  be  less  than  or  equal  to  the 
magnetosheath  density  and  to  decrease  with 
Increasing  Inward  penetration  of  the  boundary 
layer.  Thus  the  most  consistent  Interpretation 
of  the  results  places  Tordo  Uno  In  the  poleward 
portion  of  the  cleft. 

The  'high-energy'  portion  of  the  electron 
spectrum  observed  by  Tordo  but  not  by  Hawlceye  In 
Figure  6 can  be  attributed  to  an  Internal 
magnetospherlc  source,  l.e.,  drifting  energetic 
electrons  from  Che  nlghtslde.  Recent  papers  by 
McDlamild  et  al.  (1976)  and  Haerendel  and 
Paschsiann  [19731  have  shown  Chat  trapped 
energetic  electrons  are  observed  concurrently 
with  entry  layer  and  cleft  plasma.  Eastman  et 
al . 11976]  find  that  electron  spectra  In  the 

boundary  layer  develop  high-energy  (few  keV) 
calla  with  Increasing  Inward  penetration  of  the 
boundary  layer  and  have  associated  these  to 
leakage  of  magnetospherlc  plasma  Into  the 
boundary  layer.  In  the  present  case,  Tordo  Uno 
sampled  only  Inside  the  loss  cone,  and  thus 
nothing  can  be  said  concerning  the  pitch  angle 
distribution.  Spectrally,  Che  electrons  are 
s.lmllar  CO  low-altitude  polar  satellite 
measurements  In  the  same  region  [Wlnnlngham. 
19721  . 

The  Initial  energy-dependent  Increase  In 
electron  flux  with  altitude  (<250  km)  la  due  to 
the  Interaction  of  the  primary  beam  [Banks  et 
al . , 1974,  Figure  16}  with  the  atmosphere. 


Above  this  altitude,  spatial  and/or  temporal 
changes  predominate.  In  summary,  Tordo  Uno  and 
the  associated  barium  Jet  were  launched  Into  the 
poleward  portion  of  the  afternoon  cleft  which 
contained  unenergized  magnetosheachllke 
electrons.  With  reference  to  the  work  by 
Jeffries  et  al.  [1975],  Figure  3 shows  the  path 
of  the  barium  plasma  leaving  the  vicinity  of 
Tordo  Uno  to  be  essentially  directly  antisolar 
across  the  polar  cap  Into  the  nlghtslde  auroral 
zone. 

By  comparison,  we  know  the  position  of  Tordo 
Dos  vls-i-vla  the  cleft  equatorward  boundary 
relatively  well.  Tordo  Dos  was  launched 
equatorward  of  the  cleft,  crossed  Its 
equatorward  boundary,  and  encountered  two 
inverted  V's  In  the  cleft  prior  to  barium 
Injection.  In  assuming  the  cleft  boundary  to  be 
zonal  (l.e.,  parallel  to  magnetic  latitude)  over 
the  short  longitude  span  traversed  and  to  be 
fixed  for  the  duration  of  the  Tordo  Dos  flight, 
barium  injection  occurred  approximately  1.9® 
Into  the  cleft.  The  latter  assumption  Is 
supported  by  concurrent  results  from  the 
sounders. 

The  electron  precipitation  observed  In  Tordo 
Dos  was  dramatically  different  from  Tordo  Uno. 
'Monoenergetlc'  peaks  were  observed  (Figure  10b) 
In  the  electron  spectra,  and  the  energy  of  these 
peaks  varied  systematically,  l.e.,  the  inverted 
V's,  Such  behavior  has  been  attributed  to 
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Fig.  8.  Electron  energy-time  spectrogram  for  Tordo  Dos. 


acceleration  by  an  electric  field  parallel  to 
the  earth's  magnetic  field  [Burch  ct  al . . 1976, 
and  references  therein].  To  Investigate  this 
possibility,  the  numerical  model  of  Evans  [1974] 
was  applied  to  the  data.  It  should  be  noted 
that  the  experimental  data  are  Inadequate  In 
pitch  angle  and  spectral  resolution  to  prove  or 
or  disprove  uniquely  any  theory;  only 
consistency  can  be  shown. 

With  the  aforementioned  reservations  in  mind 
the  experimental  spectra  are  easily  derivable 
from  Evans'  model.  An  example  of  a spectrum 
from  the  second  Inverted  V Is  given  in  figure 
11.  A least  squares  Maxwellian  fit  Is  made  to 
the  high-energy  tall  of  the  observed  fluxes,  and 
the  temperature  of  the  assumed  Maxwellian  plasma 
at  the  top  of  the  accelerator  Is  determined  as 
indicated  on  the  figure.  The  probable  error 
from  the  least  squares  fit  Is  also  indicated. 
The  Initial  Maxwellian  is  allowed  to  fall 
through  a potential  drop  V^,  producing  a jetted 
beam  at  the  bottom  of  the  accelerator.  The  beam 
then  produces  backscattered  and  secondary 
electrons  upon  Interaction  with  the  atmosphere, 
those  upgolng  electrons  with  energies  less  than 
eV  arc  reflected  by  the  parallel  electric  field 
and  would  be  observed  by  the  Instrument  as  low 
energy  downgoing  particles.  The  Maxwellian  fit 
together  with  a choice  of  V fixes  the  density  n 
of  the  assumed  Maxwellian  plasma.  The 
accelerating  potential  Is  then  varied,  producing 
the  solid  curve  at  low  energies  In  Figure  11, 
until  a 'best  fit'  to  the  data  Is  obtained. 
Thus  for  this  model  the  only  free  parameters 


which  we  vary  are  V and  the  B Ratio,  the  ratio 
of  magnetic  field  strengths  at  the  accelerator 
exit  and  at  the  top  of  the  atmosphere. 

For  the  backscatter  calculations  we  used  the 
electron  impact  ionization  cross  section  for 
atomic  oxygen  from  Ranks  et  al.  [1974]  (see 
also  Opal  et  al.  11971]).  Use  of  the  cross 
section  for  N-  would  increase  the  low-energy 
fluxes  by  about  12%  at  about  100  eV, 

Figure  12  presents  the  results  of  these 
calculations  for  the  second  inverted  V,  All 
points  except  the  first  and  last  are  derived  or 
'fitted'  quantities.  This  first  point  (denoted 
by  a circle)  corresponds  to  the  ' nonaccelarated' 
spectrum  (Figure  lOa)  encountered  just  prior  to 
the  second  Inverted  V (see  Figure  12).  The  last 
point  (also  denoted  by  a circle)  is  from  a 
temperature  fit  only  at  that  point.  It  can  be 
seen  that  even  though  the  potential  varied 
d rama t leal ly , the  density  (assuming  source 
isotropy)  and  the  temperature  were  constant 
within  one  standard  deviation.  In  addition,  the 
'nonaccelerated'  density  and  temperature  lie 
within  the  same  bounds. 

Note  that  the  nonaccelerated  values  are  local 
quantities  and  the  derived  quantities  refer  to 
the  input  of  the  'linear  accelerator,'  as 
described  above.  Thus  the  inverted  V has  as  its 
source  population  a Maxwellian  with  • density  of 
0.R7  +0.25  cm  and  a temperature  of  117  + 24 
eV. 

If  all  the  asstimptlons  in  Fvani'  model  are 
valid  in  this  case,  the  approximate  position  of 
the  exit  of  the  accelerator  can  be  determined 
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Fig.  9.  Schematic  of  geographic  location  of  observed  auroral  arcs,  Inverted  V events, 
barium  release,  and  returns  observed  by  the  lonosondes.  All  positions  have  been  re- 
duced to  the  160-km  altitude  level. 


TORDO  DOS 

IIJAN  I97S  T,>002S:02UT 
T- 319.98  10  9l9.ISltc 

SEOGSAPHIC  'coordinates 
ALTITUDE- 907km  UAT  ■ 72*  LONG -120.9* 

PITCH  ANGLE -49* 


Fig.  10a.  Electron  spectrum  obtained  prior  to 
entry  Into  second  Inverted  V observed  by  Tordo 
Dos. 


from  the  B Ratio.  The  derived  value  of  12.5 
gives  an  approximate  lower  limit  of  IR^  above 
Tordo  or  2 R_  geocentric.  This  value  is 
consistent  wltn  the  lower  limit  value  obtained 
by  assuming  that  the  colllmatlon  cone  had 
expanded  Just  enough  to  Include  the  measured 
pitch  angle. 

Our  results  are  In  basic  agreement  with  those 
of  Burch  et  al.  (1976)  except  for  one  feature. 
Burch  et  al.  pointed  out  that  In  all  the  cleft 
Inverted  V's  they  observed  the  characteristic 
energy  (temperature)  went  up  as  the  acceleration 
potential  Increased.  They  Interpreted  this  as  a 
heating  as  well  as  an  acceleration  of  the  source 
population.  The  results  as  shown  In  Figure  12 
Indicate  that  no  statistically  Important  heating 
occurred  for  this  second  Inverted  V event  on 
Tordo  Dos.  However,  some  remarks  on  the  nature 
of  the  high  energy  tall  are  In  order.  We  found 
that  with  few  exceptions  the  more  high  energy 
data  points  we  Included  In  the  least  squares 
routine,  the  hotter  the  temperature  of  the 
plasma.  Thus  the  high  energy  tall  In  the 
Inverted  V event  Is  non-Maxwelllan.  Since  the 
potential  is  chosen  primarily  to  find  a fit  to 
the  lower  energy  portion  of  the  spectrum  and 
since  this  depends  mostly  on  the  fluxes  near  the 
peak  of  the  beam,  we  used  only  two  or  three  data 
points  near  the  beam  maximum  to  determine  the 
temperature  for  the  results  shown  In  Figure  12 
(for  example.  In  Figure  11,  fluxes  at  energies 
of  634.4,  1033,  and  1655  eV).  In  addition,  data 
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Fig.  10b.  Representative  electron  spectra  ob- 
tained In  second  Inverted  V observed  by  Tordo 
Dos. 

points  at  the  lower  Intensity  levels  are 
experimentally  less  well  determined. 

To  test  this  apparent  disagreement  with  Burch 
et  al.,  we  ran  new  'fits’  using  the  criterion 
that  we  would  use  all  data  in  the  tall  with  flux 
values  greater  than  10'*  cm”  s”  sr”  eV”  . By 
using  this  criterion,  only  8 of  the  23  cases  are 
changed.  (In  one  case,  at  t-  0030:30,  the 
temperature  went  from  54  to  127  eV,  but  we  could 


Fig.  11.  Experimental  and  modeled  electron 
spectrum  for  eecond  Inverted  V in  Tordo  Doe 
data.  Model  spectrum  generated  by  using  the 
numerical  model  of  Evans  [1974]. 

Li  


no  longer  find  a reasonable  'fit'  to  the 
spectrum.)  The  temperatures  of  these  eight  cases 
did  cone  out  higher  than  before,  necessitating 
reduced  values  for  In  some  cases,  with  n 

general  trend  of  temperature  and  V 
proportionality.  Linear  regression  analysis  In 
this  case  yields  a temperature  of  0.06  V + 102 

eV,  but  the  correlation  Is  not  high  (correlation 
coefficient  la  0.35). 

It  therefore  appears  that  our  analysis  is 
consistent  with  no  apparent  heating  of  the 
plasma  at  energies  near  the  beam  maximum,  but 
some  heating  at  higher  energies  proportional  to 
V cannot  be  ruled  out.  It  is  equally  possible 
that  the  non-Maxwell Ian  nature  of  the  tall  Is  a 
feature  of  the  original  plasma  (see  Figure  10a). 

Burch  et  al.  [1976]  reported  the  cleft 
Inverted  V's  to  be  In  a region  of  antlsunward 
convection.  No  convection  measurements  were 
available  at  the  position  of  the  Tordo  Dos 
Inverted  V's.  However,  the  barlian  streak  was 
released  approximately  7 km  poleward  of  the 
second  inverted  V and  drifted  antlsunward  at 
approximately  constant  magnetic  latitude 
[Jeffries  et  al..  1975,  Figure  3]. 

In  conclusion,  Tordo  Dos  crossed  the 
equatorward  boundary  of  the  cleft,  traversed 
approximately  1.9°  of  its  extent,  and 
encountered  t%ro  inverted  V's.  The  inverted  V's 


Fig.  12.  Density  and  temperature  derived  from 
Evans'  [1974]  stodel  for  source  spectrum  of 
second  Inverted  V observed  by  Tordo  Dos  and 
the  necessary  potential  drop.  See  text. 
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can  be  well  modeled  by  the  model  of  Evans  (197A] 
by  utilizing  a constant  source  temperature  and 
density  and  a varying  potential.  A 'best  guess' 
places  the  accelerator  at  a minimum  geocentric 
distance  of  2 Rg.  The  barium  Injected  just 
subsequent  to  the  second  Inverted  V drifted 
essentially  parallel  to  the  cleft  (i.e.,  at 
approximately  constant  magnetic  latitude)  and 
antlsunward  as  opposed  to  the  perpendicular  and 
antisunward  motion  of  Tordo  Uno.  The  cloud 
motion  and  its  dynamics  will  be  detailed  in  a 

• subsequent  paper. 
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Time-dependenl  behavior  of  tubes  of  F layer  plasma  is  computed  for  tubes  carried  around  several  flow 
paths  in  the  polar  region.  The  flow  paths  are  those  proposed  previously  by  Knudsen  (1974).  Ionization 
sources  include  direct  and  scattered  solar  photons  and  measured  fluxes  of  precipitating  energetic  elec- 
trons. Computed  electron  concentrations  are  compared  with  measured  concentrations  from  topside 
sounder  data  obtained  in  the  same  satellite  pass  as  the  energetic  electron  flux  data.  The  following 
conclusions  are  drawn:  The  proposed  convection  pattern  produces  a tongue  of  F layer  plasma  over  the 
polar  cap  with  electron  concentrations  consistent  with  the  measured  concentrations.  Had  the  F layer  been 
assumed  to  be  nonconvecting  over  the  polar  cap,  the  computed  concentrations  would  have  been  a factor 
of  10  too  small.  Rapid  convection  of  plasma  across  the  cleft  prevents  significant  increase  in  /V.F,  at  the 
cleft.  Rapid  convection  must  exist  through  the  intense  nightside  auroral  energetic  particle  precipitation 
zones,  or  a compensating  electron  loss  mechanism  must  develop  to  prevent  the  buildup  of  an  ionization 
ridge.  A mid-latitude  trough  is  formed  by  the  proposed  convection  pattern  with  only  normal  F layer 
recombination  processes  operating.  The  low  concentration  in  the  trough  is  maintained  by  scattered  solar 
EU  V photons.  The  trough  predicted  by  the  numerical  model  is  not  observed  in  the  topside  sounder  data. 


Introduction 

In  a previous  paper,  Knudsen  [1974]  presented  a model  for 
the  convection  field  of  the  high-latitude  F layer  and  evaluated 
semiquantitatively  the  expected  time-dependent  behavior  of  a 
tube  of  F layer  plasma  carried  around  the  polar  regions  by  the 
field.  The  proposed  field  appeared  to  explain  many  of  the  high- 
latitude  F layer  anomalies.  The  purpose  of  the  present  paper  is 
to  present  the  initial  results  of  a more  detailed  numerical  study 
of  the  behavior  of  the  high-latitude  F layer.  Tubes  of  ioniza- 
tion are  subjected  to  time-dependent  ionization  rates  from 
both  solar  photons  and  precipitating  energetic  electrons  as  the 
tubes  follow  the  convection  held  proposed  in  the  previous 
paper,  and  the  time-dependent  response  of  the  plasma  within 
the  tube  is  computed. 

Numerical  Model 

The  time-dependent  behavior  of  the  plasma  within  a mag- 
netic flux  tube  was  computed  with  a numerical  code  developed 
by  Schunk  and  Walker  [1973].  The  code  solves  the  coupled 
momentum  and  continuity  equations  for  NO*,  Oi*,  and  O* 
ions  in  the  E and  F regions.  Minor  changes  were  made  in  the 
code  to  permit  the  addition  of  ionization  rates  resulting  from 
energetic  electron  fluxes  and  also  from  solar  EUV  scattered 
into  the  nightside  of  the  ionosphere.  The  code  was  also  modi- 
fied to  permit  the  base  of  the  flux  tube  to  change  in  latitude 
and  longitude  with  elapsed  time  in  a prescribed  manner. 

The  code  does  not  solve  the  energy  balance  equation,  and 
consequently,  it  is  necessary  to  specify  T,  and  T,.  For  this 
study  we  have  set  both  T,  and  T|  equal  to  the  neutral  gas 
temperature  T„.  Heating  of  the  thermal  electron  gas  by  fluxes 
of  soft  electrons  with  energies  of  a few  tens  of  electron  volts  is 
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thus  neglected,  and  the  enhancement  of  electron  concentra- 
tions at  altitudes  above  approximately  400  km  in  response  to 
electron  heating  will  not  be  reproduced  in  the  model. 

The  variations  of  the  neutral  constituents  Ni,  O,,  and  O with 
latitude  and  local  time  were  derived  from  the  model  atmo- 
sphere of  Jacchia  [1971]  for  winter  (solar  declination,  -23”) 
and  medium  solar  activity  = 130  X 10""  W m"*  Hz"'). 
The  NO  density  was  the  same  as  that  used  by  Schunk  and 
Walker  [J973]. 

The  numerical  solution  was  obtained  over  the  altitude  range 
I20-S00  km.  The  upper  altitude  limit  was  set  at  S(X)  km  to 
ensure  proper  convergence  of  the  solution  for  reasonably  sized 
time  steps.  The  fluxes  of  NO*  and  O,*  were  set  equal  to  zero  at 
the  upper  boundary.  To  simulate  the  loss  of  O*  through 
operation  of  the  polar  wind,  the  upward  velocity  of  O*  ions  at 
the  SOO-km  boundary  was  assumed  to  be  constant  at  10  m s"'. 
At  the  lower  boundary  the  0«*  and  NO*  concentrations  were 
arbitrarily  set  to  l(P  cm"’.  Also  at  this  boundary  the  O* 
density  was  assumed  to  be  in  chemical  equilibrium  and  was  set 
equal  to 

' ' *.[N.]  + it,[0,] 

where  [O*  ],  [Nt],  and  [Ot]  are  the  appropriate  densities  and 
/^O*)  the  ion  production  rate  of  O*;  ilci  and  k,  are  reaction 
rates  defined  by  Schunk  and  Walker  [1973].  The  concentra- 
tions of  NO*  and  Oi*  near  the  lower  boundary  are  governed 
primarily  by  chemical  reaction  rates,  so  that  their  concentra- 
tions a short  distance  above  the  lower  boundary  are  insensitive 
to  the  arbitrarily  assigned  values. 

The  gyrofrequency  of  the  ions  is  comparable  to  or  less  than 
their  collision  frequency  at  and  below  an  altitude  of  approxi- 
mately 180  km.  The  ions  do  not  'follow*  the  flux  tubes  below 
this  altitude.  We  may  expect  therefore  that  the  model  results 
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Pig.  1.  Production  rates  of  O’,  Nj*.  and  O,*  by  solar  photons  at 
three  SZA’s  computed  by  the  numerical  model.  Those  at  109°  SZA 
represent  background  rates  from  scattered  photons. 


below  180  km  will  be  in  error  at  the  £ and  Ft  altitudes  near 
boundaries  separating  regions  of  different  ion  production 
rates. 

The  vertical  drift  which  may  be  imparted  to  the  plasma  by 
the  electric  field  causing  the  convection  of  the  plasma  and  also 
by  neutral  winds  has  been  neglected  in  this  first  model.  Vertical 
drift  was  felt  to  be  a second-order  effect  which  could  be  studied 
after  first-order  effects  were  clarified.  Some  discussion  of  the 
expected  drift  is  given  by  Knudsen  [1974]. 

The  displacement  of  the  magnetic  pole  from  the  geographic 
pole  imparts  a universal  time  dependence  into  the  high-lati- 
tude F layer  concentration,  which  is  to  be  the  subject  of  a 
future  study.  For  the  present  study  we  assume  that  the  mag- 
netic and  geographic  axes  are  collinear. 

The  ionization  sources  included  in  the  model  are  solar  EUV 
photons,  both  direct  and  scattered,  and  precipitating  electrons. 
The  method  of  calculating  the  direct  photo-ionization  rates  for 
Nj.  Oj,  and  O is  described  by  Schunk  and  Walker  [1973].  We 
are  interested  in  the  high-latitude  F layer  during  winter  in  the 
present  study  and  have  calculated  the  Chapman  grazing  in- 
cidence function  using  -23°  for  the  solar  declination  angle.  To 
simulate  the  ionization  rates  from  scattered  EUV,  the  ioniza- 
tion rates  from  direct  photons  were  added  to  background 
ionization  rates  produced  by  scattered  photons.  The  back- 
ground ionization  rates  were  computed  with  the  same  code 
that  was  used  for  the  direct  ionization  rates  with  the  following 
changes.  The  solar  zenith  angle  was  set  to  zero,  and  the  photon 
fluxes  in  the  1 1 spectrum  intervals  suggested  by  Hinlereggeret 
at.  [1965]  were  set  to  zero  except  for  the  wavelength  intervals 
1027-91 1 A,  630-460  A,  and  370-280  A.  The  fluxes  in  these 
intervals  were  set  at  I X 10",  I X 10',  and  I X 10*  photons  cm"* 
s '.  respectively,  and  correspond  to  H Lyman  0,  He  I 485  A, 
and  He  Lyman  a [Chen  and  Harris,  1971].  The  total  photo- 
ionization rates  from  both  the  direct  and  the  scattered  photons 
for  three  solar  zenith  angles  (SZA's)  are  given  in  Figure  I. 
Those  for  a SZA  of  109°  are  effectively  the  background  rates 
with  no  contribution  from  the  direct  rates. 

Ionization  rates  from  precipitating  electrons  were  computed 
by  using  energy  distribution  functions  measured  by  the  soft 
particle  spectrometer  (SPS)  on  the  Isis  2 satellite.  The  energy- 
time  spectrograms  over  the  northern  polar  cap  on  December 


15  (day  349,  orbit  3278),  1971,  are  shown  in  Figures  2 and  3. 
The  satellite  crossed  the  northern  polar  cap  from  approxi- 
mately 1 200  M LT  (magnetic  local  time)  to  2300  M LT  (Figures 
2 and  3).  Spectra  at  times  labeled  l-l  I were  selected  as  being 
representative  and  were  used  in  computing  the  ionization  rates 
produced  by  the  electrons.  The  locations  of  the  satellite  at  the 
times  at  which  the  spectra  were  recorded  are  shown  in  Figure 
7.  The  electron  precipitation  was  assumed  to  be  uniform  in 
local  time  for  approximately  3 hours  on  either  side  of  the 
satellite  track  in  both  the  dayside  cleft  region  and  the  nightside 
auroral  region.  Auroral  oval  emission  data  in  5577  A and  3914 
A recorded  on  the  Isis  2 satellite  [Lw  and  Anger,  1973]  indicate 
that  on  the  nightside  of  the  oval  the  precipitation  was,  in  fact, 
more  intense  in  the  premidnight  sector  than  in  the  post- 
midnight sector. 

The  data  from  the  satellite  pass  shown  represent  a quiet 
period.  The  AE  index  was  below  100  y for  28  hours  prior  to 
the  pass  except  for  two  short  excursions.  Sixteen  hours  prior  to 
the  pass  the  AE  index  went  as  high  as  300  y during  a 2-hour 
period,  and  approximately  3 hours  before  the  pass  it  rose  to 
approximately  150  y for  a 2-hour  period. 

The  proton  number  and  energy  fluxes  over  the  polar  region 
for  the  pass  shown  were  2 orders  of  magnitude  lower  than  the 
corresponding  electron  fluxes  and  have  been  neglected  as  an 
important  source  of  ionization. 

Electron  spectra  for  several  pitch  angles  in  the  downgoing 
hemisphere  at  each  of  the  1 1 regions  enumerated  in  Figures  2 
and  3 have  been  analyzed  and  used  to  derive  the  spectra 
illustrated  in  Figures  4 and  5.  The  precipitating  flux  was  rea- 
sonably uniform  with  pilch  angle,  and  the  spectra  shown  are 
considered  representative  of  the  spectra  observed  at  each  re- 
gion. The  SPS  on  Isis  2 measures  electron  flux  down  to  5 eV. 
For  the  purpose  of  computing  ionization  rates  from  these 
fluxes  we  have  arbitrarily  extended  the  measured  10-eV  fluxes 
to  I eV  at  a constant  level.  Since  the  numerical  code  used  for 
computing  thermal  plasma  behavior  does  not  balance  energy 
and  since  electrons  with  less  than  10-eV  energy  produce  no 
significant  ionization,  this  extrapolation  has  no  significant  ef- 
fect on  the  numerical  results. 

The  ion  production  rates  from  the  electron  spectra  pre- 
sented above  were  computed  with  an  early  version  of  a com- 
puter code  for  the  interaction  of  energetic  electrons  with  the 
atmosphere  developed  by  Banks  et  al.  [1974].  Production  rates 
for  regions  2 and  9,  representative  of  the  cleft  and  nightside 
auroral  regions,  are  shown  in  Figure  6. 

The  flow  paths  around  which  the  plasma  tubes  were  carried 
in  the  numerical  simulation  are  those  presented  in  an  earlier 
paper  by  Knudsen  [1974]  and  are  reproduced  as  Figure  7, 
Results  were  obtained  only  for  paths  labeled  Il-VI.  The  flow 
paths  define  the  motion  of  plasma  tubes  as  seen  by  an  observer 
in  a nonrotating  frame  of  reference  looking  down  on  the  north 
polar  cap.  That  is,  the  motion  includes  corotation  so  that  the 
exposure  of  the  plasma  tubes  to  solar  EUV  photons  is  properly 
simulated.  The  tubes  of  plasma  were  convected  across  the 
dayside  cleft  at  a velocity  of  I km/s  and  on  across  the  cap  at  a 
velocity  of  0.5  km/s.  They  were  again  convected  across  precip- 
itation region  7 (Figure  3)  at  1 km/s.  Equatorward  of  the  limit 
of  closed  field  lines  the  tubes  were  convected  between  succes- 
sive dots  in  I hour.  Analysis  of  energetic  electrons  measured 
on  Isis  2 at  greater  than  20  keV,  greater  than  40  keV,  and 
greater  than  200  keV  indicates  that  the  field  lines  were  closed 
equatorward  of  69° (16  min,  40  s.  Figure  3)  and  that  no 
trapped  distributions  existed  poleward  of  70° in  the  night 


I ncrgy-imic  spectrogram  lor  electrons  from  the  Isis  2 soft  partjcie  spcclTomeler.  December  UHl.  local  noon  portion  I ocations  at 
vvhtch  electron  spectra  v^ere  analv/ed  for  use  in  the  model  are  labeled 
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Fig.  4.  Representative  electron  particle  flux  as  a function  of  energy. 

Numbers  refer  to  locations  labeled  in  Figure  2. 

scctor(l6  min,  20  s,  Figure3)(J.  R.  Burrows,  private  commu- 
nication, 1976).  The  limit  of  closed  field  lines  in  the  night 
sector  of  Figure  7 is  consistent  with  this  observed  boundary. 
The  tubes  were  convected  slowly  through  precipitation  regions 
8-10  as  required  by  the  model,  and  as  we  shall  see,  these 
precipitation  regions  produced  considerable  ionization  in  the 
slowly  convecting  tubes. 

The  flow  pattern  of  Figure  7 in  the  dayside  cleft  region  and 
over  the  polar  cap  is  reasonably  consistent  with  that  derived 
from  vector  ion  velocity  data  by  Heelis  et  al.  (1976].  Flow 
across  the  cleft  and  into  the  polar  cap  is  observed  between 
approximately  0900  and  1500  MLT  as  has  been  assumed 
herein,  and  the  general  flow  pattern  of  Heelis  et  al.  presented 
with  corotation  removed  in  their  Figure  8 would,  with  corota- 
tion  restored,  look  similar  to  that  of  our  Figure  7.  The  Harang 


Fig  5.  Representative  electron  particle  flux  as  a function  of  energy. 
Numbers  refer  to  locations  labeled  in  Figure  3. 


Fig.  6.  Production  rates  of  O*,  N,*,  and  O,*  ions  by  the  electron 
particle  fluxes  at  locations  2 and  9 of  Figures  2 and  3. 


discontinuity  and  stagnation  points  must  exist  in  any  reason- 
able model  at  approximately  the  locations  indicated  in  Figure 
7,  and  the  flow  must  approach  corotation  in  the  vicinity  of 
SS’A.  The  major  discrepancy  between  experimental  measure- 
ments and  the  'time  average'  representation  presented  in  Fig- 
ure 7 appears  to  be  in  the  magnitude  of  the  velocity  equa- 
torward  but  in  the  near  vicinity  of  the  limit  of  the  closed  field 
lines.  The  measured  velocity  is  of  the  order  of  I km/s,  whereas 
the  model  velocities  with  corotation  removed  are  of  the  order 
of  0. 1 km/s.  The  measured  velocity  does  approach  corotation 
velocity  somewhere  between  TO^A  and  50*A  [Heelis  et  al., 
1976).  The  velocities  with  which  the  plasma  tubes  were  trans- 
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LIMIT  or  CLOSED  FIELD  LINES 

e STAGNATION  POINT 

PLASMA  PAUSE 

Fig.  7.  Flow  paths  around  which  time-dependent  plasma  response 
was  computed.  Locations  of  electron  spectra  are  indicated  by  the 
numbered  dots.  Numbers  correspond  to  those  in  Figures  2 and  3. 
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Fig.  8.  Contours  of  N^F,  derived  from  the  numerical  analysis.  The 
short-dashed  portions  of  contours  are  extrapolated  contours. 

ported  across  the  cleft  and  polar  cap  are  reasonably  consistent 
with  observation. 

The  finite  time  steps  in  the  numerical  analysis  were  less  than 
or  equal  to  S min.  In  regions  where  the  ion  production  rate  was 
changing  rapidly,  the  time  step  was  decreased  appropriately. 

Results 

Contours  derived  from  the  numerical  results  are 

shown  in  map  view  in  Figure  8.  The  maximum  electron  con- 
centration in  each  vertical  profile  around  the  several  flow  paths 
was  plotted  and  subsequently  contoured.  The  short-dashed 
portions  of  the  contours  equatorward  of  the  lowest  latitude 
path  on  the  nightside  of  the  earth  are  extrapolated.  Numerical 
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Fig.  9.  Synoptic  N„F,  contours  for  the  Antarctic  polar  region 
(after  Sato  and  Rourke  1 196.1]  and  reproduced  by  Thomas  and  Andrews 
1 1969]).  Coordinates  are  geographic  coordinatn. 


Fig.  10.  Comparison  of  numerically  computed  electron  concentra- 
tion with  that  measured  by  the  topside  sounder  experiment  on  Isis  2 in 
vertical  profile. 

computations  for  flow  paths  equatorward  of  the  plasmapause 
were  not  performed  for  this  study.  Figure  9 shows  experimen- 
tal data  included  for  later  comparison. 

A local  noon-midnight  vertical  profile  comparison  of  the 
numerically  computed  electron  concentration  with  the  ob- 
served electron  concentration  is  shown  in  Figure  10.  The  ob- 
served data  are  from  the  topside  sounder  experiment  on  Isis  2 
for  the  same  pass  as  that  from  which  the  energetic  electron  flux 
data  were  derived.  Comparison  of  the  model  results  with  the 
topside  sounder  data  is  meaningful  in  the  altitude  interval 
350-500  km. 

In  Figure  1 1 we  show  a comparison  of  the  computed  elec- 
tron concentration  in  vertical  profile  with  the  plasma  tube 
convecting  and  with  it  stationary  (steady  state)  at  two  loca- 
tions along  path  IV.  The  observed  electron  concentration  is 
also  indicated.  Figure  12  shows  the  computed  electron  concen- 
tration at  midnight  on  path  II.  The  plasma  within  the  con- 
vecting tube  in  this  region  of  path  II  is  in  a steady  state 
condition. 

Discussion 

The  computed  electron  concentration  shown  in  map  view  in 
Figure  8 shows  several  of  the  features  characteristic  of  the 
polar  ionosphere.  First,  a ‘tongue'  of  plasma  extends  from  the 


Fig.  II.  Comparison  of  numerically  computed  electron  concentra- 
tion profiles,  convecting  and  steady  state,  in  the  cleft  (location  I ) and 
within  the  polar  cap  (location  2)  with  measured  concentration.  Loca- 
tion I corresponds  to  77*A,  IIS4  MLT.  I(X)*  SZA,  and  particle 
spectrum  2.  Location  2 corresponds  to  84°A,  2I(X)  MLT,  117*  SZA, 
and  particle  spectrum  5. 
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local  noonside  of  the  earth  across  the  polar  cap  to  the  night- 
side.  Manifestations  of  this  tongue  in  time-averaged  f,Ft  data 
[Fike,  1971  ] and  in  Alouettc  I topside  sounder  data  at  3S0-km 
altitude  [Nishida,  I967|  were  pre.scnted  in  the  prr/ious  paper 
[Knudsen,  19741.  In  Figure  9 we  have  reproduced  synoptic 
data  for  the  Antarctic  polar  region  in  geographic  coordi- 
nates [5aro  and  Rourke.  1963;  reproduced  by  Thomas  and 
Andrews,  1969).  A tongue  is  clearly  visible  with  low  concentra- 
tion regions  on  either  side.  The  tongue  and  low  concentration 
regions  on  either  side  are  reproduced  in  the  numerical  results. 
Not  present  in  Figure  9 are  the  ridge  of  ionization  produced  by 
auroral  electron  precipitation  and  the  deep  trough  present  in 
the  numerical  results.  These  discrepancies  will  be  addressed 
hereafter.  The  region  of  low  concentration  on  the  local  evening 
side  of  the  tongue  in  the  numerical  model  results  from  the 
presence  of  a stagnation  point  outside  the  region  of  particle 
precipitation  (Figure  7)  and  direct  solar  photon  production. 
Near  the  stagnation  point  in  local  winter  the  tubes  of  plasma 
move  sufficiently  slowly  fur  the  concentration  to  decay  to  low 
values.  The  region  of  low  concentration  on  the  morning  side  of 
the  tongue  in  the  numerical  model  results  from  lack  of  signifi- 
cant exposure  of  the  plasma  tubes  to  the  sun  and  to  the  more 
intense  nightside  auroral  electron  precipitation  at  latitudes 
below  70° A in  the  present  model. 

The  direct  solar  photons  are  able  to  produce  a significant  F 
layer  for  solar  zenith  angles  from  0°  up  to  approximately  104° 
(compare  Figure  I ).  On  the  midnight  side  of  the  104°  SZA  of 
Figure  8.  only  Ihe  scattered  photons  make  a contribution  to 
the  ion  production. 

In  the  vme-averaged  data  presented  by  Nishida  [1967]  a 
maximum  existed  in  the  evening  sector  at  about  60°A.  The 
cause  of  the  maximum  is  uncertain,  although  Nishida  [1967] 
assumed  that  it  was  a manifestation  of  the  mid-latitude  eve- 
ning increase  [Evans.  1965].  The  data  used  fur  the  averaging 
included  data  from  August  29  to  November  10,  so  that  the 
stagnation  point  would  have  been  illuminated  by  direct  solar 
photons  at  nearly  all  universal  times.  The  maximum  is  not 
evident  in  the  data  presented  by  Pike  [1971]  or  by  Saio  and 
Rourke  [1964]. 

No  significant  increase  In  was  evident  in  Ihe  numerical 
results  in  the  cleft  region  (Figures  8 and  10)  As  has  been 
pointed  out  previously  [Knudsen.  1974),  a large  increase  can- 
not be  expected  if  the  plasma  is  convected  across  the  cleft  at  I 
km/s.  For  the  present  example  the  maximum  increase  to  he 
expected  is 

SN„F,  ptr  =•  2 X 10^  cm  ’ s ' x 120  s 2.4  X 10*  cm  • 

The  value  Is  only  1041'  of  N^Fi  prior  to  entry  of  the  plasma  into 
Ihe  cleft.  Ffad  the  electron  concentration  been  in  a steady  slate 
in  the  cleft  region.  Ihe  value  of  N„F,  would  have  been  a factor 
of  4 larger  in  the  cleft  than  just  equatorward  of  the  cleft 
(Figure  II).  The  observed  electron  concentration  at  400-km 
altitude  does  show  a small  increase  in  concentration  al  Hie 
cleft.  We  a.ssume  this  to  be  primarily  a result  of  the  heating  of 
the  electron  gas  as  the  lube  convects  across  the  cleft  The 
increase  in  observed  concentration  in  the  cleft  is  small  at  400 
km  but  increases  with  altitude.  This  behavior  indicates  an 
increase  in  scale  height  and  not  an  increase  in  The 

increase  would  be  independent  of  altitude  if  the  electron  and 
'on  temperature  were  held  constant,  as  is  evidencrd  in  Ihe 
lumerical  results  of  Figure  1 1 . As  has  been  pointed  out  pre- 
viously. Ihe  numerical  model  did  not  include  energy  balance, 
and  the  electron  and  ion  temperatures  were  arbitrarily  set 
equal  to  Ihe  neutral  temperature. 


Fig.  12.  Computed  electron  concentration  protile  at  2400  MLT  on 
path  II.  The  convccling  and  steady  state  solutions  are  the  same.  The 
observed  electron  concentration  Is  also  shown. 


The  computed  electron  concentration  compares  favorably 
with  the  topside  sounder  electron  concentrations  in  the  alti- 
tude interval  in  which  the  two  sets  of  data  overlap.  The  10  X 
10*  cm"’  contours  of  both  sets  of  data  coincide  rather  well 
from  75°A  on  the  dayside  to  75° .\  on  the  nightside  of  the 
model.  The  prominent  variations  in  concentration  exhibited 
by  the  topside  sounder  data,  especially  at  higher  altitude,  are 
produced  by  increases  in  Ihe  electron  and  ion  temperatures. 
This  is  evidenced,  as  was  previously  mentioned,  by  the  de- 
crease in  electron  concentration  variation  as  the  altitude  de- 
creases. The  model  results  are  not  expected  to  reproduce  these 
variations  because  the  electron  and  ion  temperatures  were  held 
constant.  The  increase  in  concentration  of  the  model  results  at 
67°A  in  the  night  sector  resulted  from  the  increase  in  ion 
production  rate  produced  by  energetic  electrons  (regions 
8-10).  The  concentration  trough  at  ~57®A  results  from  a lack 
of  production  and  a long  decay  time.  These  latter  features  do 
not  occur  in  Ihe  experimental  data. 

The  nearly  constant  value  of  N„Ft  over  the  polar  cap 
evidently  results  from  rapid  convection  of  plasma  from  the 
dayside  of  Ihe  cleft.  The  concentration  profiles  that  would 
occur  without  convection  (steady  slate)  are  shown  in  Figure  1 1 
and  are  compared  with  the  profiles  resulting  from  convection. 
The  nonconvecting  profile  at  liKation  2 is  an  upper  limit 
profile  computed  from  the  production  rates  of  ions  for  the 
upper  limit  particle  precipitation  flux  labeled  5 in  Figures  2 
and  4 The  profile  would  be  identical  to  the  steady  stale  profile 
in  Figure  12  were  the  particle  flux  zero. 

Because  Ihe  polar  cap  F layer  statistically  has  a large  con- 
centration and  decays  acrou  Ihe  cap  from  noon  to  midnight  in 
a manner  consistent  with  convection  [Nuhido.  I%7;  Knudsen. 
1974).  It  IS  reasonable  to  infer  that  F layer  plasma  must  be 
convected  rather  steadily  across  the  cleft  near  the  local  noon- 
time sector  Were  convection  of  ionospheric  plasma  across  the 
cleft  and  on  across  the  polar  cap  to  cease  for  a period  of 
several  hours,  the  polar  cap  F layer  would  decay  by  as  much  as 
an  order  of  magnitude. 

The  ridge  of  high  ion  concentration  occurring  in  the  numer- 
ical results  equatorward  of  70°  on  the  nightside  of  the  earth 
(Figure  8)  is  produced  by  assumed  slow  convection  of  plasma 
through  the  intense  particle  precipitation  zones  represented  by 
spectra  8 and  9 of  Figure  3.  This  ridge  of  ionization  is  not 
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present  in  the  observed  electron  concentration  (Figure  10). 
Statistical  studies  also  show  little  or  no  enhancement  in  the 
nightside  auroral  /one  [sec  Knudsen.  1974],  Suppression  of 
such  a ridge  could  be  accomplished  by  moving  the  plasma 
through  these  /ones  rapidly.  Rapid  convection  would  not  per- 
mit buildup  of  ioni/ation  m a tube  as  it  crossed  the  zone(s)  of 
high  particle  precipitation.  The  buildup  may  also  be  sup- 
pressed through  an  increased  loss  rate  within  the  auroral  zone 
[5('Aun4:  el  ai.  I976|. 

A trough  region  of  low  ion  concentration  was  produced  in 
the  numerical  results  equatorward  of  the  ridge  of  high  ion 
concentration  from  approximately  1800  MLT  to  0600  MLT. 
The  trough  resulted  from  the  slow  convection  of  plasma  near 
the  evening  stagnation  point  at  which  location  the  solar  zenith 
angle  was  too  large  for  direct  solar  photons  to  maintain  the  F, 
layer  ion  concentration  at  a high  value  against  the  normal  ion 
loss  proces.ses  and  polar  wind  escape.  A ‘steady  state'  level  of 
ion  concentration  of  flow  line  II  around  the  nightside  of  the 
polar  region  was  maintained  by  the  scattered  solar  photon  flux 
level  assumed  to  be  present  in  the  model.  The  trough  dis- 
appears on  the  morning  side  as  the  plasma  tube  emerges  into 
the  direct  solar  photon  flux.  The  ionization  level  in  the  trough 
would  be  least  near  midnight  and  greater  toward  the  evening 
and  morning  sectors  in  a more  realistic  model  in  which  the 
Kattered  solar  photon  flux  decreased  with  increasing  solar 
zenith  angle. 

The  deep  trough  evident  in  the  numerical  results  (Figures  8 
and  10)  did  not  appear  in  the  topside  sounder  data.  The 
probable  explanation  is  that  variation  of  the  convection  elec- 
tric held  and  hence  flow  paths  with  time  prevents  the  develop- 
ment of  the  deep  trough.  The  flow  paths  must  be  steady  for 
periods  of  the  order  of  24  hours  for  a deep  trough  to  develop  in 
the  manner  suggested  by  Knudsen  [1974]  and  demonstrated  in 
the  present  numerical  results.  As  was  pointed  out  earlier,  the 
magnetic  AE  index  was  small  and  steady  except  for  two  short 
excursions  for  28  hours  prior  to  the  orbit  selected  for  analysis. 
Hence  the  conditions  for  development  of  a trough  by  slow 
convection  would  seem  to  have  been  fulfliled  as  well  as  one 
could  expect.  Resolution  of  this  discrepancy  must  await  future 
studies. 

Summary  and  Conclusions 

The  time-dependent  behavior  of  the  F layer  plasma  within  a 
magnetic  flux  tube  is  computed  for  the  tube  curried  around 
several  flow  paths  in  the  polar  regions.  The  flow  paths  are 
those  proposed  in  a previous  paper  by  Knudsen  [1974].  As  the 
tube  traverses  each  path,  it  is  subjected  to  time-dependent 
ionization  rates  from  direct  and  scattered  solar  photons  and 
precipitating  energetic  electrons.  Ionization  rates  for  the  ener- 
getic electrons  are  computed  from  energy  spectra  observed  by 
the  Isis  2 soft  particle  spectrometer.  The  numerical  results  are 
presented  in  the  form  of  a map  view  of  fV«Ft  contours,  elec- 
tron concentration  in  vertical  section  over  the  magnetic  pole 
from  noon  to  midnight,  and  several  vertical  profiles  of  electron 
concentration  at  several  locations  for  both  convecting  and 
nonconvecting  flux  tubes.  The  computed  electron  concentra- 
tions arc  compared  with  observed  electron  concentrations 
from  topside  sounder  data  obtained  in  the  same  satellite  pass 
as  the  energetic  electron  flux  data. 

The  proposed  convection  field  produced  a tongue  of  F layer 
plasma  extending  from  the  dayside  of  the  cleft  over  the  polar 
cap  with  concentrations  consistent  with  the  observed  concen- 
trations. Had  the  F layer  been  assumed  to  be  stationary  (non- 
convectmg)  within  the  polar  cap.  the  computed  concentrations 


would  have  been  a factor  of  10  too  small.  Concentration  ‘lows’ 
on  either  side  of  the  tongue  are  consistent  with  some  pre- 
viously observed  synoptic  data.  No  significant  increase  in 
NmF,  occurred  at  the  cleft  in  the  model  results.  This  behavior 
is  consistent  with  the  observed  behavior.  In  diffusive  equilib- 
rium an  increase  of  approximately  4 in  concentration  would 
have  existed. 

The  presence  of  the  polar  cap  F layer  with  peak  electron 
concentrations  typically  in  excess  of  10“  cm  ' implies,  when  it 
is  interpreted  in  the  light  of  the  present  study  and  a previous 
study  [Knudsen.  1974],  that  F layer  plasma  is  rather  steadily 
convected  across  the  cleft  in  the  local  noontime  sector  and  on 
across  the  polar  cap. 

A ridge  of  ionization  not  present  in  the  observed  ionosphere 
was  predicted  by  the  numerical  simulation  in  the  nightside 
auroral  zone.  Several  zones  of  particle  precipitation  were  mea- 
sured by  the  Isis  2 satellite,  and  the  plasma  tubes  were  trans- 
ported at  high  velocity  (I  km/s)  only  through  the  highest 
latitude  zone.  Evidently,  the  plasma  tubes  are  transported 
rapidly  across  zones  of  high  precipitation,  and/or  some  addi- 
tional loss  mechanism  is  operative  in  the  precipitation  zones. 

A trough  of  low  electron  concentration  was  predicted  by  the 
numerical  model  which  extended  from  approximately 
1800-2400  MLT  to  approximately  0600  MLT.  The  electron 
concentration  in  the  trough  was  constant  with  local  time  and 
maintained  by  the  assumed  constant  flux  of  solar  EUV  pho- 
tons scattered  into  the  night  ionosphere.  No  trough  was  pre- 
sent in  the  Isis  2 topside  sounder  data.  A satisfactory  ex- 
planation for  the  latter  discrepancy  must  await  future  studies. 
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The  Development  of  Auroral  and  Geomagnetic  Substorm  Activity  After  a 
Southward  Turning  of  the  Interplanetary  Magnetic  Field 
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A comprehensive  study  of  growth  phase  and  substorm  activity  following  a period  of  magnetic  calm  has 
been  conducted  through  a network  of  all-sky  camera  stations,  auroral  zone  magnetic  observatories,  and 
panicle  detectors  aboard  the  Isis  1 satellite.  We  have  carefully  documented  the  observations  with  the 
following  results.  The  preexpansive  phase  arc  which  extended  at  least  from  17  to  05  MLT  was  responsible 
for  an  energy  input  rate  of  =<3  X 10“  ergs/s  before  breakup.  An  equatorward  drift  of  this  arc  of  6 
km/min.  observed  only  in  the  evening  sector,  remained  until  after  the  expansive  phase,  when  its  motion 
stopped  abruptly  at  the  time  of  the  maximum  poleward  displacement  of  the  arcs.  Electrons  responsible  for 
the  prebreakup  arc  had  energies  of  at  1-5  keV.  Protons  of  a4-keV  energy  were  measured  equatorward  of 
the  electron  arc.  During  the  expansive  phase,  symmetrically  traveling  disturbances  were  observed  propa- 
gating eastward  in  the  evening  sector  and  westward  in  the  morning  sector.  The  propagation  stopped  for 
1-2  min  at  the  time  of  maximum  expansion  and  then  continued,  thus  suggesting  a momentary  variation  m 
the  rate  of  convection.  Equivalent  currents  consistent  with  observed  magnetic  perturbations  represented 
approximately  the  same  DPZ  (twin  vortex)  pattern  before  the  expansive  phase  as  during  it;  however, 
although  the  magnitude  of  the  currents  was  greater  during  the  expansive  phase,  the  dominant  feature 
during  this  phase  was  an  intense  westward  auroral  electrojet.  The  camera  observations  of  diffuse  cloudlike 
aurora  showed  an  injection  of  iv40-keV  electrons  during  the  expansive  phase  along  the  auroral  oval 
between  midnight  and  0400  corrected  geomagnetic  time.  Movement  of  the  cloud  indicated  an  eastward 
gradient  drift  of  the  electron  population. 


Introduction 

On  the  basis  of  data  from  an  extensive  body  of  all-sky 
camera  documentation  and  polar  magnetograms,  Akasofu 
[1964]  outlined  the  substorm  sequence.  Although  his  initial 
description  has  been  modified  slightly  [Akasofu,  1968:  Mont- 
brtand,  1971],  it  remains  essentially  correct  in  ordering  the 
observed  geophysical  phenomena  into  a consistent  picture. 
Since  that  time  an  enormous  effort  has  gone  into  further 
documenting  substorm-related  phenomena  and  into  estab- 
lishing the  cause  of  the  onset  of  the  substorm  expansive  phase 
(see  review  by  Rostoker  \\912\).  Whereas  in  Akasofu's  descrip- 
tion the  substorm  sequence  is  initiated  by  the  sudden  breakup 
of  an  auroral  arc  accompanied  by  a steep  decrease  of  the 
horizontal  component  of  the  magnetic  field,  it  was  claimed  by 
McPherron  [1970]  that  the  complete  substorm  sequence  in- 
cludes a disturbance  period  which  he  denoted  the  substorm 
growth  phase  According  to  his  statement,  intervals  of  mag- 
netic calm  may  be  followed  by  significant  deviations  of  the 
horizontal  component  of  the  magnetic  field  prior  to  the  start  of 
the  expansion  phase  of  a magnetospheric  substorm. 

Since  McPherron's  suggestion  that  a growth  phase  exists, 
several  authors  [liiima  and  Sagata.  1972;  McPherron  et  ai. 
1973;  Kokuhun  and  liiima.  1975]  (and  others)  have  contributed 
to  the  study  of  the  phenomenon  by  giving  more  detailed  de- 
scriptions of  the  signatures  which  are  observed  to  be  character- 
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istic  of  that  situation.  The  authors  appear  to  agree  in  the 
opinion  that  the  onset  of  an  isolated  substorm  is  preceded  by  a 
growth  phase,  which  is  closely  related  to  the  north-to-south 
change  of  the  interplanetary  magnetic  field  (IMF),  following  it 
with  a delay  of  10-20  min.  The  signatures  of  this  growth  phase 
observed  on  the  ground  are  the  following:  a gradual  decrease 
in  H,  the  horizontal  magnetogram,  at  auroral  zone  stations 
before  its  sharp  drop;  a gradual  decrease  in  H at  low  latitudes 
(especially  in  the  evening-midnight  sector);  the  growth  of  a 
polar  equivalent  current  of  the  twin  vortex  mode  (especially  in 
the  polar  cap):  and  the  equatorward  motion  of  auroras  [Koku- 
bun  and  lijima,  1975], 

In  opposition  to  the  above  mentioned  view,  Akasofu  et  al. 
[1973]  have  suggested  that  the  changes  in  the  magnetic  field 
characteristic  of  the  growth  phase  have  limited  significance, 
since  these  changes  do  not  themselves  result  in  the  onset  of  the 
expansive  phase.  They  suggested  instead  that  the  mechanism 
which  triggers  the  expansive  phase  is  an  internal  one,  largely 
independent  of  external  factors  related  to  convection  such  as 
the  southward  turning  of  the  interplanetary  magnetic  held.  In 
continuation  of  this  suggestion,  Akasofu  [1975]  has  demon- 
strated that  the  southward  turning  of  the  B,  component  and 
the  subsequent  chain  of  processes,  proposed  by  McPherron  et 
al.  [1973],  do  not  have  any  significance  in  causing  the  ex- 
pansive phase.  Substorms  are  frequently  being  triggered  in 
direct  relation  to  such  a chain  of  processes,  in  which  the 
‘growth  phase'  of  the  substorm  may  then  appear  as  part  of  the 
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TABLE  I.  Stations  Used  in  This  Study 

Geographic  Corrected 

Geomagnetic 


Station 

Symbol 

Latitude, 

deg 

Longitude, 

deg 

Latitude, 

•N 

Observation 

Alert 

AT 

82.5 

295.5 

86 

M 

Thule 

TH 

77.5 

290.8 

86 

A,  M 

Resolute  Bay 

RB 

74,7 

265.2 

84 

M 

Godhavn 

GO 

69.2 

306.5 

77.5 

A,M 

Ny  Aalesund 

78.9 

11.9 

75.5 

A 

Sdndre  Strdmfjord 

SS 

67.0 

309.2 

75,5 

A 

Baker  Lake 

BL 

64.3 

264.0 

75 

M 

Heiss  Island 

BT 

80.6 

58.5 

74.5 

M 

Fort  Churchill 

FC 

!8.8 

265.9 

71 

A 

Narssarssuaq 

NAR 

61.2 

314.6 

68 

A.M 

Great  Whale  River 

GW 

55.3 

282.2 

68 

A.M 

Dixon  Island 

Dl 

73.5 

80.4 

67.5 

M 

T romsd 

TR 

69.7 

18.9 

66.5 

A.  M 

Reykjavik 

RY 

64.2 

338.3 

66 

M 

Tixi  Bay 

Tl 

71.6 

129,0 

65 

M 

Murmansk 

MM 

69.0 

33.1 

64.5 

S.  M 

Kiruna 

Kl 

67.8 

20.4 

64.5 

M 

College 

CO 

64.9 

212.2 

64.5 

M 

Meanook 

ME 

54.6 

246.6 

64 

M 

Wellen 

WE 

66.2 

190.2 

60 

M 

Lerwick 

LE 

60.1 

358.8 

59 

M 

Bangui 

BA 

4.4 

18.6 

5 

M 

M'Bour 

MB 

14.6 

343.0 

21 

M 

San  Juan 

SJ 

18.1 

293.8 

30 

M 

Fredericksburg 

FR 

38.2 

282.6 

50 

M 

T ucson 

TU 

32.2 

249.2 

40 

M 

Honolulu 

HO 

21.3 

202.0 

21 

M 

complete  substorm  picture.  However,  they  may  as  well  be 
released  at  a later  stage,  after  the  end  of  the  B,  negative 
interval,  as  long  as  the  area  of  the  polar  cap  is  greater  than  a 
certain  minimum  value,  namely,  the  area  after  a prolonged 
period  of  large  positive  B,.  Once  the  polar  cap  has  achieved 
this  minimum  value,  no  substorm  can  be  released  until  a new 
southward  turning  of  the  IMF  has  initiated  erosion  of  field 
lines  on  the  dayside  of  the  magnetosphere  with  the  accom- 
panying features  of  increasing  convection,  observable  from 
ground-based  observatories. 

No  matter  whether  the  growth  phase  is  an  integral  part  of 
the  substorm  phenomenon,  as  proposed  by  MePherron,  or 
more  likely  an  independent  B,  negative,  increased  convection 
phenomenon  in  accordance  with  the  statements  of  Akasofu, 
there  is  still  a need  for  further  detailed  study  of  the  develop- 
ment of  the  geomagnetic  and  auroral  activity  observed  after 
the  southward  turning  of  the  IMF. 

On  February  25,  1969,  the  Isis  I satellite  happened  to  pass 
over  the  premidnight  sector  of  the  auroral  oval  during  a 
geophysical  situation  which  appears  to  have  been  identical 
with  the  substorm  growth  phase  defined  by  MePherron.  The 
event  followed  a period  of  low  substorm  activity,  characterized 
by  values  of  /tf  less  than  200  y during  the  34  hours  before  the 
substorm  and  less  than  60  y during  the  last  6 hours.  By 
combining  the  satellite  particle  measurements  with  ground- 
based  observations  we  have  been  able  to  report  in  the  follow- 
ing a careful  documentation  of  the  development  of  the  auroral 
and  magnetic  activity  through  the  growth  phase  and  the  sub- 
sequent substorm  onset  and  expansion  phase  as  well  as  of  the 
particle  precipitation  preceding  the  substorm  onset.  Our  ob- 
servations give  a detailed  picture  of  the  gradual  development 
of  auroral  and  geomagnetic  activity  which  follows  the  south- 
ward turning  of  the  interplanetary  magnetic  field  and  which,  in 


this  particular  event,  is  completed  by  the  onset  of  a substorm 
with  remarkable  poleward  expansion. 

Data  Sources 

In  our  study  we  made  use  of  ground-based  observations 
made  during  the  event  at  the  stations  listed  in  Table  I , together 
with  particle  measurements  from  the  Isis  I satellite  and  mag- 
netic field  measurements  from  the  satellites  Explorer  35  and 
Heos  I.  The  ground-based  observations  are  auroral  all-sky 
photographs  (A  in  column  6 of  Table  1),  spectrograms  (S  in 
column  6),  and  geomagnetic  records  (M  in  column  6). 

The  interplanetary  magnetic  field  was  measured  by  the  satel- 
lites Explorer  35  in  the  evening  quadrant  ({.X,  Y,  Z)»e  = (-8. 
-t-63,  -61)  Be)  and  Heos  I in  the  forenoon  quadrant  ((Af,  F, 
2)(jhm  * (II,  -31,  14)  Be).  Both  satellites  appear  to  have  been 
outside  the  shock  front.  None  of  the  records  was  complete; 
thus  detailed  results  from  Heos  I are  missing  during  00(X>^300 
(JT,  but  hourly  averages  do  exist  for  this  interval  [A(/r|g,  1975]. 
Similarities  of  gross  features  of  the  B,  variation  at  the  two 
satellites  indicate  a time  lag  between  them  of  about  30  min. 
From  the  position  of  the  satellites  it  is  estimated  then  that  the 
B,  variations  at  the  front  of  the  magnetopause  occur  15-17 
min  earlier  than  those  at  Heos  I.  Hence  the  B,  variations 
observed  at  Heos  I have  been  referenced  to  the  front  of  the 
magnetosphere  by  a shift  of  15  min  toward  earlier  hours.  In 
Figure  I they  have  been  shown  together  with  magnetic  records 
from  the  stations  Tromsd,  Reykjavik,  Narssarssuaq,  Great 
Whale  River,  and  Fort  Churchill,  which  were  all  situated  in 
the  18-04  magnetic  time  sector  of  the  auroral  zone.  In  the 
figure  the  stations  are  ordered  from  east  to  west.  Each  arrow 
represents  2(X)  y.  The  magnetic  H records  from  a number  of 
mid-  and  low-latitude  stations  are  shown  in  Figure  2.  The 
stations  are  ordered  from  east  to  west,  with  the  most  easterly 
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Fig.  I.  Normal  magnetograms  from  the  auroral  zone  stations 
Iromsi).  Reykjavik,  Narssarssuaq.  Great  Whale  River,  and  Fort 
Churchill,  February  24-25,  1969.  Each  arrow  represents  a variation  of 
2(X)  y.  The  dots  indicate  station  geomagnetic  midnight.  Top  curve:  B, 
(vertical  GSM  coordinate)  as  measured  al  Heos  I,  referenced  to  the 
front  of  the  magnetosphere  by  a shift  of  1 5 min  toward  earlier  hours. 


one  (Bangui)  at  the  top  of  the  figure.  Each  arrow  in  the  figure 
represents  20  >. 

At  the  time  of  the  substorm  the  nightside  of  the  auroral  zone 
was  situated  between  western  Siberia  in  the  east  and  North- 
west Territories  in  the  west,  with  the  magnetic  midnight  meri- 
dian crossing  the  Iceland-Greenland  area.  Conditions  were  not 
ideal  for  auroral  observations.  The  moon  was  in  its  third 
quarter;  the  Soviet  cameras,  Reykjavik,  Kiruna,  and  several 
stations  in  Greenland  were  not  in  operation  because  of  bad 
weather,  and  most  of  the  remaining  stations  had  some  cloud 
cover  early  in  the  evening.  Nevertheless,  eight  well-distributed 
cameras  were  available  for  a study  of  the  development  of  the 
auroral  display  in  the  course  of  the  night  (Table  I and  Figure 
3).  At  each  of  the  stations  the  sky  was  photographed  with  one 
exposure  per  minute  except  at  Tromsd,  where  three  exposures 
per  minute  were  taken.  The  photographs  have  been  used  as  a 
basis  for  the  description  of  the  distribution  of  the  aurora;  a 
selected  series  has  been  transformed  to  synoptic  maps  in  geo- 
graphical coordinates,  shown  here  as  Figures  4<j  and  4b.  In 
projecting  the  auroral  forms  on  the  maps  a mean  height  of 
1 10-120  km  has  been  assumed  for  the  lower  border  of  the  arcs 
and  bands. 


Magnetic  Observations 

The  IMF  had  been  directed  northward  for  several  hours 
when  the  decrease  of  B,  began  between  2200  and  2300  UT. 
Because  of  the  oscillatory  character  of  the  change  it  is  difficult 
to  determine  the  exact  hour  when  B,  went  negative;  after 
2345-2350,  B,  was  definitely  negative.  The  first  departure  from 
several  hours  of  magnetic  quiet  was  observed  at  approximately 
the  same  time,  i.e.,  at  2345-(X)l5  at  the  auroral  zone  stations 
situated  near  and  shortly  after  magnetic  midnight  (Figures  I 
and  5).  The  departure  is  seen  as  a gradual  decrease  of  H,  but  D 
and  Z are  also  affected.  In  the  evening  sector.  Great  Whale 
and  Fort  Churchill  were  little  influenced  by  the  disturbance 
before  01(X)  UT  (Figures  I and  6).  At  the  low-latitude  stations 
situated  in  the  evening  sector  a gradual  decrease  of  H set  in 
between  (XXX)  and  OOiO  UT  and  continued  until  the  expansive 
phase  onset  at  0130.  The  onset  is  recognized  as  a sudden 
decrease  of  H in  Figure  1 . In  Figure  2 it  appears  as  the  onset  of 
a positive  bay  in  the  midnight  sector  and  of  a negative  bay  of 
comparable  magnitude  in  the  evening-afternoon  sector.  The 
transition  takes  place  near  the  longitude  of  San  Juan.  At  San 
Juan  a flat  depression  is  observed  from  shortly  after  midnight 
to  03(X)  UT,  with  only  a small  negative  perturbation  to  be  seen 
in  connection  with  the  expansion  phase  of  the  substorm.  The 
gradual  decrease  of  H at  mid-  and  low-latitude  stations  as  well 
as  the  flat  depression  observed  at  San  Juan  may  be  caused  by 
an  increase  of  a cross-tail  or  ring  current  system,  which  be- 
comes detectable  in  the  records  about  j hour  after  the  change 
of  sign  of 

The  magnetic  signatures  of  this  event  are  in  agreement  with 
the  description  given  by  MePherron  [1970].  Following  a de- 
crease of  B,  to  negative  values  a growth  phase  of  about  U-hour 
duration  is  observed,  characterized  by  a gradual  decrease  of  H 
at  low  latitudes  and  by  a contemporary  change  of  the  magnetic 
elements  at  auroral  latitudes  in  the  direction  in  which  they  are 
subsequently  changed  by  the  substorm.  The  growth  phase  is 
superseded  by  the  expansion  phase  initiated  by  the  sharp  onset 
at  about  0130. 

Before  the  end  of  the  recovery  phase  of  this  substorm  a new 
onset  is  noticed  at  0210-0215  at  the  low-latitude  stations, 
especially  at  Bangui  and  M’Bour.  In  the  auroral  zone  the 
magnetic  bay  is  most  clearly  recorded  at  Narssarssuaq;  at 
Tromsd  the  onset  is  recognized  as  an  abrupt  change  of  the 
slope  in  H.  The  magnetic  activity  from  this  double  substorm 
dies  away  by  about  0400  UT.  After  2 hours  calm  a new 
moderate  substorm  is  observed  in  the  midnight  sector  (Fort 
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Fig.  2.  Horizontal  variation  at  low-  and  mid-latitude  stations, 
February  24-25,  1969.  Length  of  arrows  is  20  T.  The  dots  indicate 
station  midnight;  the  open  circle  at  Honolulu  indicates  noon. 
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Development  of  presubstorm  arc,  February  2S,  1969.  The  times  indicated  are  UT.  Circles  represent  the  fields  of 
view  of  camera  stations  (A)  of  Table  I. 


Fig  ia  Development  of  auroral  display,  February  25,  l%9.  Growth  phase  and  poleward  expansion 
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Churchill,  Figure  I ).  We  shall  comment  further  on  the  onset  of 
this  substorm  in  a later  section. 

We  have  used  the  ground  magnetometer  data  in  the  con- 
struction of  equivalent  current  vectors  to  represent  the  current 
distribution  before  and  after  breakup  at  01 30  UT.  This  enables 
a comparison  of  the  current  system  during  the  growth  phase 
with  that  during  the  time  of  maximum  expansion.  We  have 
taken  the  time  of  the  satellite  crossing  of  the  arc  (0 1 1 5 UT ) as  a 
time  during  which  we  would  consider  the  growth  phase  to  be 
well  in  progress.  The  equivalent  current  vectors  which  could 
produce  the  magnetic  perturbations  at  a number  of  polar  cap 
and  auroral  zone  stations  at  the  time  of  the  crossing  have  been 
drawn  in  a corrected  geomagnetic  latitude/time  diagram  in 
Figure  7.  Also  shown  are  the  equivalent  current  vectors  de- 
rived from  the  same  stations  during  the  time  of  maximum 
disturbance  of  the  event.  The  maximum  magnetic  disturbance 
was  measured  during  the  20  min  immediately  following  the 
poleward  expansion:  the  situation  at  0150  UT  has  been  chosen 
for  construction  of  the  current  vectors.  The  quiet  level  for  both 
graphs  has  been  defined  as  the  level  which  immediately  pre- 


ceded the  onset  of  the  event,  i.e.,  during  2200-2400  UT.  The 
equivalent  current  system  during  the  growth  phase  is  similar  to 
the  system  corresponding  to  fi,  < -I  y in  the  winter  months 
(DPZ  system)  as  constructed  by  Friis-Christensen  and  Wil- 
hjelm  [1975].  In  this  system,  which  is  believed  to  depict  the 
magnetic  field  produced  by  field-aligned  current  sheets  in  the 
auroral  oval  and  by  a Hall  current  maintained  by  the  electric 
field  between  the  sheets,  the  disturbance  vectors  over  the  polar 
cap  are  of  the  same  order  of  magnitude  as  those  in  the  post- 
midnight auroral  oval.  In  Friis-Christensen  and  Wilhjelm's 
paper  (their  Figure  5)  the  equivalent  vectors  are  approximately 
2 times  the  magnitude  of  those  found  here,  in  accordance  with 
the  fact  that  their  parameter  6,  (the  2-hour  average  of  the 
actual  and  preceding  hour)  is  numerically  greater  in  their  data 
than  in  ours. 

Following  a further  decrease  of  B,  the  equivalent  currents 
during  the  maximum  disturbance  at  01 50  UT  are  more  intense 
than  those  at  01 15  by  a factor  of  2,  thus  being  comparable  to 
the  average  currents  of  Figure  5 of  Friis-Christensen  and  H'lV- 
hjelm  (1975).  Superimposed  on  the  DPZ  system  we  find  at 


Fig.  4*.  Development  of  auroral  display.  February  25.  19W.  Recovery  and  second  (minor)  substorm 
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Fig  5.  Rapid  run  magnetograms  showing  the  horizontal  com- 
ponent at  Narssarssuaq  and  Reykjavik.  February  24-25,  1%9.  The 
gradual  depression  prior  to  abrupt  onset  of  the  magnetic  substorm  is 
clearly  displayed. 


0150  UT  an  intense  westward  auroral  electrojel.  Thus  al- 
though the  current  systems  during  growth  {B,  negative)  and 
expansive  phase  are  identical  in  most  respects,  they  differ  by 
the  presence  during  the  expansive  phase  of  the  intense  west- 
ward auroral  electrojet.  A similar  statement  has  been  given 
earlier  by  lijima  and  Sagaia  [1968].  The  dominant  role  of  the 
westward  electrojet  after  breakup  is  in  agreement  with  the 
substorm  model  of  McPherron  et  al.  [1973],  according  to 
which  the  breakup  is  associated  with  the  sudden  increase  of  the 
westward  electrojet. 

A URORAL  Observations 

The  course  of  the  auroral  substorm  near  magnetic  midnight 
can  be  summarized  by  reference  to  Figure  8.  The  figure  pre- 
sent.s,  in  negative,  a sequence  of  all-sky  photographs  from 
Narssarssuaq.  The  breakup  occurred  at  0128  UT,  corrected 
geomagnetic  midnight  at  the  station.  The  poleward  expansion 
ended  at  about  0137;  after  this  time  the  oval  part  of  the  bulge 
was  gradually  replaced  by  a diffuse  veil,  leaving  the  discrete 
forms  in  the  polar  front  of  the  bulge  only.  The  auroras  re- 
turned from  the  north,  and  by  0213  a new  expansion  was  in 
progress  to  the  west  of  the  station,  giving  rise  to  increased 
intensity  in  the  withdrawing  bands.  After  the  second  expan- 
sion, auroral  arcs  remained  at  about  71°  corrected  geomag- 
netic latitude  until  a new  substorm  on.set  occurred  near  0600 
UT.  It  follows  from  the  magnetic  and  auroral  observations 
presented  above  that  the  period  studied  may  be  divided  into 
two  time  intervals  with  rather  different  types  of  activity,  the 
prehreakup  interval  following  the  calm  period  at  about  0045 
and  the  postbreakup  (expansion  and  early  recovery)  interval 
which  started  at  about  0130  UT.  In  the  following  sections  we 
present  and  discuss  in  more  detail  ground-based  observations 
made  during  the  event  at  the  stations  listed  in  Table  I,  for  the 
two  intervals  separately  Observations  of  drift  of  auroral  arcs 
and  surges  which  were  made  within  both  intervals,  spanning 
the  lime  of  breakup,  are  presented  in  a special  section,  as  are 
also  the  particle  measurements  from  the  Isis  I satellite. 

Prebreakup  Interval 

Buildup  of  Aclivilv 

The  first  observation  of  auroral  light  before  the  breakup  was 
made  near  the  midnight  meridian  at  Narssarssuaq  station. 


Although  it  is  difficult  to  be  certain  because  of  the  presence  of 
moonlight  and  a few  clouds,  a faint  arc  stretching  from  the 
northeastern  horizon  up  to  about  10°  elevation  is  possibly  seen 
in  the  frames  taken  at  2353,  2358,  and  later  at  0006-(XX)8  UT. 
From  the  time  of  0013-0014  an  auroral  band  was  definitely 
present  and  extended  from  the  northeast  horizon  to  the  zenith. 
The  band  persisted  as  a single  homogeneous  arc  until  at  0050 
the  intensity  increased  gradually,  beginning  at  the  eastern  hori- 
zon. Folds  shaped  like  a narrow  westward  traveling  surge 
moved  from  the  horizon  toward  the  zenith  and  combined  with 
the  original  arc  at  0100  to  form  a row  of  multiple  arcs  of 
moderate  intensity.  A few  minutes  later  a new  'surge’  was  seen 
low  in  the  northeast  for  about  5 min.  The  intensity  increases 
and  the  westward  movement  of  folds  were  accompanied  by  a 
negative  bay  in  H with  an  amplitude  of  20  y in  Narssarssuaq 
and  Reykjavik  (Figure  5). 

Apparently  as  a continuation  of  the  formation  of  parallel 
arcs  over  Narssarssuaq,  the  arc  developed  to  the  west,  becom- 
ing first  detectable  within  the  Great  Whale  field  of  view  at  0100 
UT,  at  the  eastern  horizon  of  Fort  Churchill  at  0106,  and 
crossing  the  sky  to  the  northwestern  horizon  at  0109  (Figure 
3).  At  Narssarssuaq,  after  a few  minutes  of  intensity  decrease, 
a sharp  arc  developed  from  the  horizon  in  the  southwest  and 
extended  to  the  eastern  horizon.  During  the  same  time  the 
display  was  developing  toward  the  east.  It  was  observed  from 
Tromsd  through  a light  cloud  cover  at  01 16  UT  as  a homo- 
geneous arc  reaching  from  the  western  horizon  toward  the 
northeastern  horizon;  at  0117  it  stretched  from  horizon  to 
horizon.  The  arc  could  have  been  present  in  lower  intensity 
before  01 16,  so  that  the  observation  might  indicate  an  increase 
in  intensity  in  an  already  existing  aurora.  Certainly,  by  0125 
UT  the  arc  system  extended  at  least  from  a geographic  longi- 
tude of  100°W  to  45°E  (17.5-05.5  MLT)  with  sufficient  in- 
tensity for  easy  camera  observation  of  the  luminosity. 

At  about  0115  UT  a poleward  turning  of  the  arc  was  ob- 
served in  low  elevation  at  the  eastern  horizon  of  Great  Whale 
(Figure  18).  The  arc  was  first  observed  at  Great  Whale  at  0100. 
Between  this  time  and  01 15  the  arc  resembled  that  of  Figure  18 
without  the  poleward  turning  at  the  extreme  east.  The  exact 
minute  at  which  the  distortion  sets  in  is  difficult  to  determine, 
because  the  region  is  insufficiently  covered  by  the  cameras.  For 
the  same  reason  It  is  impossible  to  judge  whether  the  arc  is 
continuous,  thus  forming  either  a poleward  bulge  or  a spiral  in 
this  area,  or  discontinuous,  forming  a sun-aligned  arc  as  a 
signature  of  increased  convection  (a  short  polar  cap  band 
becomes  visible  at  approximately  the  same  time  over  Spitsber- 
gen (Figure  3)).  In  the  following  we  refer  to  this  part  of  the 
display  as  ‘the  poleward  extension.’  Further  signatures  of  the 
formation  of  the  poleward  extension  may  also  have  been  the 
temporary  increase  in  the  intensity  of  the  western  part  of  the 
arc  as  observed  from  Narssarssuaq  during  the  interval 
01 13-01 16,  as  well  as  the  onset  of  small  irregular  pulsations 
with  periods  of  several  minutes  in  the  Great  Whale  magneto- 
gram (Figure  6).  The  disturbance  set  in  between  0110  and 
0115,  probably  about  0112,  and  continued  as  a gradual  de- 
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Fig  6 Rapid  run  magnetogram  of  the  horizontal  component  al 
Great  Whale  River,  February  25.  1969.  Note  the  small-scale  per- 
turbations during  formation  and  approach  of  westward  traveling 
surge. 
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Fig.  7,  Equivalent  current  vectors  during  prebreakup  (growth ) phase  (0115  UT,  left)  and  manimum  of  expansive  phase 
(0150,  right)  of  subslorm  on  February  25, 1969.  Numbers  beside  arrows  indicate  vertical  magnetic  perturbation  in  gammas. 
Note  the  different  scale  in  the  two  diagrams. 


crease  of  100  7 as  the  westward  traveling  surge  reached  the 
station  at  about  0139. 

Following  the  formation  of  the  poleward  extension  the  in- 
tensity again  increased  from  the  east  in  the  southern  arc  (which 
at  this  time  was  the  only  one  visible)  over  Narssarssuaq.  the 
lower  border  became  irregular,  and  the  intensity  increased 
from  the  west  along  the  arc  to  meet  the  easterly  increase  near 
the  meridian.  Here  the  arc  formed  a sharp  kink  at  0125,  and 
lateral  movements  were  observed  between  0125  and  0127.  The 
intensity  increase  from  the  region  of  the  poleward  extension  to 
the  west  of  the  station  was  accompanied  by  an  obvious  though 
small  decrease  in  H (Figure  5).  The  actual  breakup  began  at 
0128-0129  very  near  Narssarssuaq  as  parallel  arcs  were 
formed  along  the  northern  border  of  the  original  arc.  The 
abrupt  onset  of  the  magnetic  perturbation  in  H occurred  at 
0128  in  the  Reykjavik  magnetograms  and  less  than  i min  later 
at  Narssarssuaq  (Figure  5). 

Interpretation  of  Observations 

The  magnetic  and  auroral  observations  may  be  summarized 
as  follows.  Within  minutes  to  tens  of  minutes  after  the  shift  of 
sign  in  fi„  precipitation  of  auroral  particles  is  observed  in  the 
auroral  oval  close  to  magnetic  midnight.  The  latitude  of  the 
oval  is  68°.  A gradual  deviation  of  the  magnetic  elements  at 
auroral  zone  stations  is  initiated  at  the  same  time,  indicating 
the  presence  of  a OPZ  (convection)  current  system,  and  signa- 
tures of  increased  cross-tail  or  ring  currents  are  seen  in  the 
records  of  //  at  low-latitude  stations.  The  observations  in- 
dicate an  increasing  convection  set  up  by  the  negative  value  of 
B,.  A minor  disturbance  (double  bay)  at  about  ObTO  is  visible 
in  stations  from  the  polar  cap  to  low  latitude  This  is  probably 
a DPI  disturbance  [Nishida.  1968],  which  may  indicate  irregu- 
larities in  the  convection  rate  The  increase  in  auroral  intensity 
at  0050  followed  by  a narrow  westward  traveling  surge  and  a 
development  toward  morning  and  evening  of  the  quiet  arc  is 
consistent  with  the  description  of  a confined  substorm  given  by 
Lui  el  al.  (1975].  In  their  example  the  magnetic  effect  of  the 


substorm  was  about  20  7.  The  decrease  of  H in  our  0050-01 10 
event  is  approximately  20  7 and  shows  the  signatures  of  a 
confined  substorm.  This  disturbance  is  immediately  followed 
by  changes  in  the  auroral  pattern  which  could  be  associated 
with  the  convection  pattern,  i.e.,  a poleward  turn  of  the  arc 
before  midnight  and  a polar  cap  aurora  in  the  morning.  The 
observations  from  the  last  10  min  before  breakup  give  the 
impression  of  increasing  convection  and  precipitation  accom- 
panied by  increasing  instability  of  the  arcs,  which  finally  result 
in  the  breakup  and  expansion. 

Expansion  and  Early  Recovery 

The  rapid  poleward  expansion,  initiated  by  the  breakup, 
was  completed  by  0137  as  the  breakup  arcs  covered  the  entire 
area  from  a few  degrees  south  of  Narssarssuaq  to  Sdndre 
StrdmQord  (Figures  8 and  9).  Simultaneously  with  this  expan- 
sion the  southernmost  arcs  moved  equatorward  with  decreas- 
ing intensity;  at  the  end  of  the  poleward  expansion,  all  discrete 
forms  except  those  at  the  highest  latitudes  dissolved,  and  an 
equatorward  broadening  of  a diffuse  veil  took  place,  which 
soon  covered  most  of  the  sky.  At  the  same  time  the  discrete 
forms  at  the  poleward  edge  of  the  display  began  a slow  equa- 
torward retreat. 

The  westward  border  of  the  region  of  poleward  extension 
between  Canada  and  Greenland  appears  to  have  remained 
relatively  stationary  during  the  expansive  phase.  However,  as 
the  expansion  proceeded,  a narrow  zonal  ‘fiord’  did  develop  in 
the  low-latitude  part  of  this  westward  extension.  The  forma- 
tion of  this  fiord  is  evident  in  Figure  Aa  in  the  0130  UT  frame 
between  60°  and  70®  west  longitude.  The  Great  Whale  photo- 
graphs show  that  it  developed  into  a small  loop  which  became 
more  and  more  elongated  toward  the  west,  eventually  opening, 
thus  forming  two  bands  a few  degrees  north  of  the  original 
band.  These  bands,  shown  in  the  0139-0157  UT  frames  of 
Figures  Aa  and  46,  reached  the  Great  Whale  meridian  al  about 
0141  but  did  not  continue  into  the  Fort  Churchill  field  of  view. 
The  development  of  the  western  loop  was  accompanied  by  a 
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I !(!  9.  Piisitioii  in  corrected  geomagnetic  latitude  of  the  auroral 
•irc'  at  Naroarssuaq.  During  the  growth  phase  the  position  is  in- 
dic.itcd  h>  the  open  dots.  The  positions  of  the  northernmost  and 
■Hiulhcrnniost  arcs  during  cspan.sion  are  shown  by  the  solid  dots.  The 
h. Itched  area  is  more  or  less  filled  by  discrete  and  dilTuse  auroras.  The 
dense  hatching  indicates  the  equatorward  border  o(  the  dilTuse  veil 

gradual  decrease  in  H at  Great  Whale  (Figure  6).  Thus  the 
loop  behaved  like  the  w-eslward  traveling  surge  as  described  by 
Akusiilu  (1968)  with  the  exception  that  the  westward  border  of 
the  poleward  expanding  bulge  did  not  move  toward  the  west 
together  with  the  surge  but  remained  stationary  on  the  evening 
side  of  the  poleward  extension  of  the  arcs.  In  a single  snapshot 
of  the  auroral  display  taken  in  the  interval  0139-0157  UT  the 
loop  would  have  been  interpreted  as  multiple  discrete  auroral 
arcs  linked  with  the  western  limit  of  the  poleward  expanding 
bulge.  .A  DMSP  satellite  photograph  showing  just  this  situa- 
tion on  November  24.  1971.  has  been  published  by  Pike  and 
H'halen  [1974].  Our  study,  based  on  all-sky  camera  data,  dem- 
onstrates how  such  multiple  arcs  may  develop  from  a loop  on 
the  western  border  of  the  poleward  expanding  bulge. 

The  propagation  toward  the  morning  sector  of  the  distur- 
bance associated  with  the  breakup  and  expansion  phase  may 
be  observed  from  Tromsd,  which  was  situated  4 hours  after 
geomagnetic  midnight.  A secjuence  of  all-sky  photographs  il- 
lustrating the  development  of  the  di.splay  at  this  station  is 
presented  m Figure  10  The  intensity  of  the  arc  immediately  to 
the  north  of  Tromsd  began  to  increase  at  0130.  Two  minutes 
later  a parallel  arc  of  lower  intensity  became  visible  on  the 
equatorward  edge  of  the  original  one  and  slowly  increased  in 
intensity,  .At  about  013t  ihe  two  arcs  were  dislocated  equa- 
lorward.  .ind  after  0139.  dilTuse  auroras  with  bandlike  struc- 
ture appe.ired  on  the  equatorward  side  of  the  arcs.  The  move- 
ment of  the  original  arc,  which  seemed  to  constitute  the 
poleward  edge  of  the  display,  as  well  as  the  movement  of  the 
boundaries  of  Ihe  dilTuse  aurora,  has  been  plotted  in  Figure  1 1. 
Ibis  figure  also  shows  the  distribution  of  auroral  light  along 
the  meridian  as  recorded  hy  the  spectral  camera  C-I80-S  at 
I oparskaya  in  the  Murmansk  area.  Auroral  data  are  not 
available  before  01.30  because  of  snow,  but  after  this  hour, 
spectrograms  of  lO-min  exposure  time  were  recorded  A com- 
parison of  the  two  plots  in  Figure  I I shows  good  agreement 
between  Ihe  stations  with  regard  to  Ihe  latitudinal  shift,  but  a 
time  lag  of  10  ± 5 min  is  observed  between  the  equatorward 
shifts.  Apparently,  Ihe  low-latitude  diffuse  area  propagated 
eastward  with  an  approximate  velocity  of  i°  'min  (or  50 
km,  min ).  On  Ihe  assumption  that  the  difru.se  cloudlike  aurora 
was  produced  by  electrons  which  had  become  trapped  on  near- 


earth field  lines  during  the  breakup  and  expansion  process  in 
the  magnetosphere,  their  energy  may  be  estimated  by  using  a 
formula  deduced  by  Sleiien  ei  at.  |I97I]  from  a paper  by 
Hamlin  el  ai  [1961].  The  formula,  which  gives  an  expression 
for  the  drift  period  of  nonrelativistic  electrons  mirroring  close 
to  the  atmosphere,  may  be  written  T(s)  = 4 x IO‘(/.£)  '. 
where  T is  the  azimuthal  drift  period  in  seconds,  E is  the 
energy  in  keV  of  the  trapped  electrons,  and  L is  the  Mcliwain 
shell  parameter.  L is  approximately  6.  and  the  stations  are 
separated  by  1 2°,  The  observed  time  lag  of  10  ± 5 min  gives  an 
estimated  particle  energy  of  £ = 37  ± 10  keV  (1^40  keV ).  If  we 
assume  that  the  electrons  became  trapped  on  near-earth  field 
lines  (£  =<  6)  no  earlier  than  at  the  breakup,  then  they  would 
have  drifted  along  this  shell  for  no  more  than  10  min  before 
arriving  at  Tromsd.  This  places  the  eastern  limit  of  the  trap- 
ping region  only  about  10°  to  the  west  of  Tromsd,  i.e.,  between 
Iceland  and  Scandinavia. 

The  observations  from  Narssarssuaq  (Figure  9)  indicate 
that  the  first  possible  incidence  of  auroral  particles  at  £ = b 
takes  place  at  0133-0134.  about  4 min  before  the  effect  at 
Tromsd.  Thus  it  seems  probable  that  the  diffuse  low-latitude 
aurora  is  formed  practically  without  time  delay  all  along  the 
auroral  oval  from  midnight  to  shortly  before  0400  geomag- 
netic time,  with  a delay  at  later  local  times  which  may  he 
explained  by  the  assumption  of  energetic  electrons  drifting  in 
the  geomagnetic  field. 

The  diffuse  cloudy  auroras  over  northern  Scandinavia  began 
to  withdraw  gradually  toward  the  auroral  zone  at  about  02(K) 
UT.  During  the  withdrawal  the  luminosity  concentrated  into  a 
few  diffuse  bands.  At  02 1 5-0220  the  southern  border  was  near 
Tromsd.  and  at  0230  the  display  was  no  longer  detectable  in 
the  all-sky  photographs.  Thus  the  poleward  retreat  of  the 
equatorward  boundary  matches  well  the  recovery  of  the  hori- 
zontal field  at  Tromsd  (Figure  1 ). 

Observation  ot- Auroral  Driet 

Equatorward  Shift  of  the  Oval 
in  the  Evening  Sector 

From  the  first  occurrence  of  auroral  light  a gradual  equa- 
lorward  drift  of  the  arc  was  observed  in  the  evening  sector. 
The  position  of  the  arc  observed  in  the  early  evening  at  f ort 
Churchill  (19  geomagnetic  time)  and  at  Great  Whale  (21  geo- 
magnetic time)  is  shown  in  Figure  12.  The  figure  shows  the 
drifts  to  be  similar  at  the  two  stations  with  a steadier  drift 
occurring  at  Fort  Churchill.  The  irregularities  at  Great  Whale 
centered  at  0122  and  0130  depict  the  passage  of  two  large 
wavelike  structures,  the  first  of  which  was  a forerunner  to  the 
small  surge  that  is  easily  detectable  in  the  photographs  of 
0121-0123  (Figure  13). 

The  equatorward  shift  of  the  arc  was  measured  to  be 
3°/hour  and  therefore  cannot  be  attributed  solely  to  the  natu- 
ral shift  of  the  position  of  the  auroral  oval  in  the  course  of  the 
evening.  The  change  in  the  latitude  of  the  Feldstein  oval  (Q  - 
3)  at  20  MLT  is  approximately  l°/hour  equatorward.  from 
which  we  can  infer  an  equatorward  drift  in  the  position  of  the 
arc  relative  to  the  oval  of  2°/hour  or  about  4 km, 'min.  This 
observation  agrees  with  results  reported  by  SlarKov  and  Eeid- 
stein  (1971)  that  auroral  forms  generally  drift  equatorward 
during  an  hour  or  so  prior  to  the  expansive  phase  of  a sub- 
storm. but  It  should  he  noted  that  in  our  obsersations  the 
equatorward  drift  was  found  only  in  the  evening  part  of  the 
oval.  Neither  at  midnight  nor  in  the  morning  sector  could  .i 
similar  drift  be  found. 
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It  is  of  parliculur  interest  to  note  that  the  eqnatoruanl  drift 
continued  through  the  growth  and  expansion  phases  but 
stopped  at  the  two  evening  stations  a few  inimitcs  before  the 
poleward  expansion  ended.  This  was  also  a few  niinutes  before 
the  narrow  westward  traveling  surge  arrixed  at  (ireat  W hale. 
If  we  make  the  tentative  assumption  that  the  eqiiatorward  drift 


of  the  arc  results  from  a westward  ionospheric  field  (see.  for 
example,  AW/cf  e/ ()/  (in'’l):  l/orer  ( 1471 1).  the  magnitude  of 
this  field,  given  by  v = (E  « B I W'.  is  4 ni\  m westward.  Our 
observations  show  that  this  tield  was  present  during  the 
growth  phase  and  most  of  the  expansion  phase  and  then 
suddenly  disappeared 
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Fig  II.  Broiidening  of  auroral  display  at  Trom.5«(  (dots)  and  Mur- 
mansk (line  segments),  based  on  all-sky  camera  (TR)  and  spectro- 
graph (MM).  The  intensity  maximum  in  the  spectrographic  record  is 
indicated  by  a narrow  double  line.  In  the  spectrogram  the  poleward 
boundary  could  not  be  determined  with  sufficient  accuracy.  After 
1)225.  observations  were  not  possible  al  T romsd  because  of  cloudiness. 


Symmetrically  Traceling  Disturbances 

In  the  morning  sector  the  equatorward  broadening  was 
initiated  by  the  passage  of  an  eastward  traveling,  small  but 
distinct  bulge  in  the  main  arcs.  The  movement  of  the  bulge  can 
be  followed  in  Figure  10.  A closer  inspection  of  the  figure 
shows  that  the  arc  was  dislocated  slightly  equatorward  on  both 
sides  of  the  bulge.  The  combined  effect  of  this  dislocation  and 
the  bulge  produces  tbe  depression  of  the  northern  arc  at 
01.1.5-0145  shown  in  Figure  1 1.  The  projection  on  the  ground 
of  the  arc  system  which  clearly  shows  the  bulge  is  presented  in 
I igure  14  More  detailed  measurements  of  the  eastward  shift 
of  the  position  of  the  bulge  are  given  in  Figure  15.  Note  that 
the  ordinate  in  the  figure  gives  the  distance  of  the  maximum 
deviation  in  the  bulge  along  the  arc  from  an  arbitrary  zero 
meridian,  measured  in  degrees  of  geomagnetic  longitude. 

At  C'lreat  Whale  a disturbance  similar  in  shape  to  but  of 
smaller  dimension  than  a westward  traveling  surge  was  .seen  to 
move  from  the  vicinity  of  the  great  poleward  bulge  between 
Canada  and  Cireenland  toward  the  west.  As  is  shown  in  the 
photographic  sequence  of  Figure  1.1.  the  surge  was  first  de- 
tected at  0120;  at  0125  it  appears  to  have  reached  the  Great 
Whale  meridian,  but  al  this  time  it  was  difficult  to  distinguish 
because  of  the  presence  of  a great  number  of  small-scale  folds 
along  the  arc.  At  01 .11-01.12  the  band  was  again  more  regular, 
with  a single  bulge  in  the  western  end.  quite  similar  to  the 
'lisliirbance  present  near  Tromsd  at  the  same  time.  The  west- 
ward movement  of  this  bulge  and  that  of  the  bulge  originally 
approaching  Great  Whale  from  the  east  (presumably  the  same 
form)  are  plotted  in  Figure  15  and  labeled  GW.  As  it  ap- 
proached the  western  horizon  at  Great  Whale,  it  became  ob- 
.servable  in  low  elevation  from  Fort  Churchill,  and  the  motion 
measurements  taken  from  the  Fort  Churchill  photographs  are 
plotted  in  the  same  figure  (labeled  FC),  A few  minutes  later  the 
f ort  Churchill  camera  slopped  operation  for  4 min.  enough 
lime  for  the  bulge  to  become  rather  poorly  defined.  But  the  few 
observations  made  at  Fort  Churchill  before  the  break  follow 
the  trend  in  the  Great  Whale  plot.  The  displacement  in  longi- 
tude results  from  the  uncertainty  in  selecting  the  same  point  on 
the  aurora  from  the  two  stations.  At  any  rale,  the  observations 
from  these  two  stations  show  that  a small  bulge  moved  west- 


ward across  the  Great  Whale  field  of  view  and  into  the  Fort 
Churchill  field  of  view  during  the  time  interval  0120-0139. 
Between  0134  and  0142  the  eastward  drift  of  a similar  bulge 
near  Trom$4  was  measured.  The  two  surges  appeared  to  have 
about  the  same  drift  speed;  in  the  time  interval  OI3S-OI38. 
when  theTromsd  measurements  were  most  accurate,  the  east- 
ward drift  speed  was  2.2°  ± 0.3°  geomagnetic  longitude  per 
minute,  while  the  westward  drift  speed  of  the  bulge  near  Great 
Whale  was  measured  to  be  2.0°  ± 0.5°  geomagnetic  longitude 
per  minute.  The  bulge  at  Tromsd  drifted  close  to  the  107° 
geomagnetic  longitude  meridian,  and  the  one  at  Great  Whale 
drifted  near  the  357°  meridian.  Thus  during  this  short  interval 
the  two  bulges  were  of  similar  shape  and  were  observed  to  drift 
symmetrically  in  relation  to  the  geomagnetic  midnight  meri- 
dian (a<50°  geomagnetic  longitude). 

The  western  irregularity  became  visible  in  the  Great  Whale 
field  of  view  a few  minutes  after  the  formation  of  the  poleward 
extension  between  Canada  and  Greenland.  The  drift  velocity 
of  this  and  accompanying  irregularities  along  the  arc  was 
relatively  constant  or  slightly  decreasing  until  0134-0135. 
when  an  increase  in  the  velocity  occurred.  Unfortunately,  data 
are  not  available  to  show  whether  the  drift  stopped  at 
0139-0140,  as  is  the  case  with  the  eastern  irregularity,  although 
the  shape  of  the  drift  curve  of  Figure  15  does  indicate  this 
possibility.  At  Tromsd  the  drift  stopped  completely  for  about 
2 min  at  CI39-0I40,  after  which  it  continued  with  the  same 
velocity  as  before.  If  the  symmetric  drift  of  the  two  irregulari- 
ties is  more  than  a coincidence,  it  can  perhaps  be  thought  of  as 
an  effect  associated  with  changes  in  the  convection  and  thus  in 
the  polar  cap  electric  field  during  the  expansion. 

Satellite  Particle  Observations 

At  0115  UT  the  Isis  I satellite  passed  over  the  arc  at  a 
position  10.3°  east  of  Great  Whale.  The  electron  spectrogram 
of  the  southbound  pass  is  presented  in  Figure  16,  where  the 
electrons  associated  with  the  visual  arc  appear  distinctly  at  an 
invariant  latitude  of  about  68°  and  a magnetic  local  time  of 
20.8  hours.  Immediately  poleward  of  this  precipitation  is  a 
restricted  region  containing  electrons  of  energies  in  the  100-  to 
250-eV  range.  Equatorward  of  the  auroral  precipitation  the 
spectrogram  shows  photoelectrons  of  energy  between  =«10  eV 
and  =>60  eV,  coming  down  the  field  lines  from  sunlit  magnet- 
ically conjugate  points. 

The  electron  diifcrential  number  spectrum  measured  near 
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Fig.  12.  Equatorward  drift  of  auroral  arc  in  the  evening  sector. 
(Top  curve)  Fort  Churchill  (19  MET).  (Bottom  curve)  Great  Whale 
River  (21  MET). 
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the  center  of  the  region  of  precipitation  is  shown  m figure  17 
This  spectrum,  laheled  d\  JE.  is  charactcrisiic  ot  nightside 
auroral  precipitation  showing  a rather  steep  negative  slope  at 
low  energies  and  the  typical  auroral  peak,  in  this  case  at  l.,1 
keV.  The  total  energy  in  the  spectrum  (in  the  range  10  eV  to 
1 1 .6  keV)  is  4.3  ergs^cm’ s sr.  .Assuming  an  isotropic  electron 
distribution,  this  energy  flur  corresponds  to  a brightness  of  20 
kR  of  light  at  5577  A [Dalgarmt  ct  al . I‘)h5|  f xamination  of 
individual  spectra  indicates  that  the  peak  in  the  electron  spec- 
trum increa.ses  and  then  decrea.scs  as  the  arc  is  crossed  Such  a 
latitudinal  morphology  has  been  referred  to  as  an  •inverted  V 
structure  by  Frank  and  Ackerson  [1971)  and  is  a common 
feature  associated  with  many  auroral  arcs 


The  enlarged  photograph  of  figure  IS  shows  the  arc  as  it 
was  observed  from  the  Great  Whale  River  Observatory  at 
0115  IJT.  The  small  white  circle  near  the  right  edge  of  the 
photograph  marks  the  point  at  which  the  geomagnetic  field 
line,  through  which  the  satellite  passed  at  the  time  the  max- 
imum energy  flux  was  detected  (01 14.5fi».  intersects  the  height 
of  maximum  auroral  luminosity.  As  is  shown  by  the  ‘energy 
flux’  curve  of  Figure  17.  most  of  the  energy  was  carried  at  ^3 
keV.  which  places  the  height  of  maximum  luminosity  at  about 
120  km  [Re«.  1969). 

The  soft  particle  spectrometer  also  provided  a measurement 
of  protons  in  the  range  of  energies  10  cV  to  1 1 6 keV.  Within 
this  energy  range,  protons  were  observed  to  precipitate  over  a 
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Fig.  14.  Maps  showing  eastward  movement  of  surge  in  the  morning  sector,  as  observed  from  Tromsd. 

region  extending  from  an  invariant  latitude  of  67.4°-66.6°.  onset  of  the  expansive  phase.  Figure  19  shows  the  electron 
The  proton  energy  flux  profile  is  plotted  in  Figure  19,  which  energy  flux  profile  measured  as  the  satellite  passed  over  the 

clearly  demonstrates  that  the  proton  precipitation  occurred  arc.  A directional  intensity  of  0.2  erg/cm’ s sr  corresponds  to  I 

equatorward  of  the  electron  precipitation  associated  with  the  kR  of  light  at  5577  A.  which  is  near  the  visibility  threshold 

visual  arc.  The  energy  associated  with  the  protons  over  this  above  the  background  light  of  the  night  sky.  As  the  figure 

range  of  invariant  latitude  averaged  3.0  X 10'*  erg/cm’ s sr.  indicates,  an  electron  flux  in  excess  of  this  value  was  detectec 

Because  of  low  count  rates  it  was  not  possible  to  measure  a for  four  consecutive  spectral  sweeps  or  for  8 s.  The  satellite 

proton  spectrum,  but  the  average  proton  energy  appeared  to  horizontal  speed  was  7 km/s,  so  that  at  the  satellite  altitude  of 
be  near  4 keV.  620  km  the  ‘arc  width’  was  56  km.  Integrating  across  the  arc. 

Since  the  satellite  crossed  the  region  of  particle  precipitation  the  rate  of  energy  input  per  unit  length  of  the  arc  was  5 X 10’ 

before  the  time  of  maximum  disturbance,  it  is  possible  to  ergs/s  cm.  The  prebreakup  arc  was  observed  by  cameras  from 

estimate  energy  input  into  the  auroral  ionosphere  before  the  west  of  Fort  Churchill  to  east  of  Tromsd  and  most  likely 


Fig  15  Drift  from  midnight  of  auroral  surges.  The  drift  seen  at  Great  Whale  and  Fort  Churchill  is  westward;  that  al 

Tromsd  is  eastward 
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Fig  17  F.leciron  difTerenlial  numher  spectrum  (r/.V  </t  ) and  en- 
ergy flux  measured  hy  the  Isis  I soft  particle  spectrometer  at  01 14  Sft 
UT  on  February  25,  1%') 


extended  much  further  in  both  directions.  A length  of  7 X 10“ 
km  approximates  this  distance,  and  with  this  assumption  the 
total  energy  input  rate  was  3.5  X 10'*  ergs^s.  The  arc  was 
visible  for  at  least  J hour  before  breakup;  therefore  an  energy 
of  7 X I0‘*  ergs  was  deposited  into  the  auroral  atmosphere  by 
the  electrons  before  breakup.  The  protons  in  the  region  equa- 
torward  of  the  arc  carried  considerably  less  energy  than  the 
electrons,  a feature  not  unusual  in  auroral  events.  The  ratio  of 
the  energy  of  the  protons  equatorward  of  the  arc  to  that  of 
electrons  within  the  arc  was  1.6  x 10  '. 

The  question  of  w hether  or  not  auroral  arcs  always  occur  on 
closed  magnetic  field  lines  continues  to  be  a matter  for  deter- 
mination We  have  attempted  such  a determination  for  the  arc 
observed  on  February  25  by  making  use  of  the  S-50  detector  of 
the  Isis  1 energetic  particle  experiment  (compliments  of  J.  R. 
Burrows,  National  Research  Council.  Ottawa.  Canada).  This 
detector  responds  to  electrons  of  energy  >40  keV  and  offers 
the  advantage  of  a large  geometric  factor  of  4 x 10  ’ cm“  sr. 
Figure  20  shows  the  S-.50  flux  plotted  against  universal  time 
and  invariant  latitude,  where  the  counts  have  been  averaged 
over  pitch  angle  in  order  to  compare  the  trapped  particle  flux 
outside  the  63°  loss  cone  with  the  downcoming  (o  < 63° ) and 
upcoming  (a  > 1 17°)  fluxes.  At  latitudes  below  about  68°  the 
trapped  flux  exceeds  the  precipitated  flux,  indicating  a closed 
field  line  configuration.  From  the  figure  it  is  not  possible  to  tell 
exactly  where  the  angular  distribution  changed,  but  it  must 
have  been  between  67.9°  and  68.4°.  Fritz  ( 1970)  would  identify 
this  location  as  .k,,.  a boundary  defined  to  be  the  lowest 
latitude  at  which  isotropy  occurs  for  electrons  of  £ > 35  keV. 
We  note  that  on  the  nightside  this  boundary  always  occurs  at  a 
latitude  lower  than  the  .V,  of  McDiarmid  and  Burrows  [1968). 
which  is  the  lowest  latitude  at  which  counts  fall  to  cosmic  ray 
background  and  which  is  interpreted  by  them  as  the  high- 
latitude  limit  of  closed  field  lines.  In  the  present  case,  X,  was 
above  latitude  71 .2°,  the  highest  latitude  tracked  on  this  pass. 
Thus  we  can  state  that  the  electron  arc  is  associated  with  or  lies 
below  the  boundary  between  open  and  closed  field  lines  and 
that  the  proton  arc  clearly  lies  on  closed  magnetic  field  lines. 


Disclssion 

(Jrowlh  Phase 

It  has  been  shown  that  ihc  subslorm  under  study  here, 
which  followed  a period  of  magnetic  calm,  appears  to  have 
been  preceded  by  an  interval  before  breakup  with  activity 
corresponding  to  the  growth  phase  defined  by  MePherron 
[1970).  The  gradual  development  of  the  activity  starts  shortly 
after  the  southward  turning  of  the  interplanetary  magnetic 
field:  the  magnetic  deflections  at  auroral  latitudes  as  well  as  at 
middle  and  low  latitudes,  the  auroral  activity,  the  westward 
magnetospheric  electric  field  deduced  from  the  auroral  drift, 
and  the  precipitation  of  auroral  particles  measured  by  the 
satellite  all  support  MePherron's  suggestion  that  the  start  of 
the  substorm  growth  phase  indicates  either  a commencement 
or  a gradual  enhancement  of  magnetospheric  convection.  The 
observation  of  symmetrically  drifting  auroral  surges  during 
growth  and  expansion  phases  as  well  as  of  auroral  patterns 
similar  to  part  of  the  convection  pattern  may  gi-e  further 
support  to  this  idea.  In  the  present  case  we  are  apparently 
concerned  with  a situation  in  which  increased  convection  initi- 
ated by  a negative  B,  is  terminated  by  several  substorms:  first  a 
very  small  confined  substorm  and  then  two  greater  substorms 
with  clear  breakup  and  poleward  expansions  which  arc  not 
separated  from  each  other  hut  which  seem  to  belong  to  the 
same  disturbance  event. 

The  very  small  substorm  occurred  at  about  OKK)  L'T  ap- 
proximately J hour  before  the  expansion  phase  of  the  major 
substorm  It  was  characterized  by  a small  westward  traveling 
surge  and  a small  magnetic  bay  of  ^2fl  y and  therefore  (its  the 
description  of  a weak  substorm  along  a contracted  oval  as 
concluded  from  ground-based  observations  \.\tonlhriand. 
1971)  and  DMSP  2 imaging  [Aka.soju.  19746),  The  southward 
turning  of  the  IMF  preceded  even  this  event  by  more  than  an 
hour,  and  the  first  departure  from  magnetic  quiet  was  ob- 
served at  auroral  /one  as  well  as  low -latitude  stations  about  I 


Fig  18  Ml-sivy  photograph  taken  at  (ireat  \V  hale  River  Observa- 
tory at  01 15  I'T  on  February  25.  IW  The  small  white  circle  desig- 
nates the  field  line  intersection  point  at  120  km 
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Fig.  19.  Electron  and  proton  energy  flux  latitudinal  profiles.  The 
proton  precipitation  lies  distinctly  equatorward  of  the  electron  precipi- 
tation associated  with  the  visual  arc.  Minute  ISofhourOI  is  labeled  as 
Oils  and  begins  at  the  00  mark  on  the  abscissa. 

hour  before  the  appearance  of  the  small  substorm.  Accord- 
ingly, we  take  the  view  that  this  substorm  may  also  have 
resulted  from  the  southward  turning  of  the  IMF  and  repre- 
sented a small  release  of  energy  before  the  much  greater  release 
at  0130. 

As  in  the  main  event  under  study,  the  0600  substorm  which 
appears  at  Fort  Churchill  and  Great  Whale  River  was  pre- 
ceded by  a short  interval  of  activity,  which  is  similar  to  the 
growth  phase  of  the  0130  substorm:  following  the  recovery  of 
substorm  activity  at  0400,  quiet  arcs  were  situated  near  71° 


corrected  geomagnetic  latitude.  Near  0525,  B,  went  negative 
(Figure  I ).  At  0525-0530  the  arcs  became  more  intense  (Figure 
21):  at  0541  an  arc  was  broadening  from  the  midnight  sector 
into  the  Fort  Churchill  field  of  view.  At  0603  a bright  aic  was 
suddenly  observed  at  about  68°:  at  0606-0610.  internal  activity 
was  seen  in  the  arc:  at  0608-0610  it  developed  toward  the  east 
across  the  Great  Whale  zenith:  and  at  0616  a moderate  pole- 
ward  expansion  began.  The  development  of  the  auroral  dis- 
play has  its  counterpart  in  the  magnetic  record  from  Fort 
Churchill  (Figure  1 ).  Also,  the  low-latitude  stations  near  mid- 
night (Fredericksburg  and  Tucson.  Figure  2)  show  features  of 
increasing  cross-tail  current  between  0540  and  substorm  onset 
at  about  0600.  Thus  although  the  observational  data  are  more 
limited  at  this  substorm,  we  see  features  very  much  like  those 
described  above,  beginning  with  the  negative  turning  of  B,  and 
ending  with  the  poleward  expansion  about  50  min  later.  The 
poleward  expansion  does  not  reach  as  high  a latitude  as  in  the 
0130  substorm:  the  product  of  B,  and  time  from  southward 
turning  to  breakup  is  also  le.ss  in  the  0600  substorm  [Arnoldy. 
1971]. 

Our  observations  do  not  determine  whether  the  expansive 
phase  trigger  mechanism  is  internal  or  external.  We  can  only 
point  to  our  observations  that  the  southward  turning  of  the 
IMF  was  followed  by  clear  signatures  of  increased  magneto- 
spheric  convection:  DPZ  equivalent  current  systems  of  gradu- 
ally (though  not  monotonically)  increasing  intensity,  signa- 
tures of  an  increasing  cross-tail  and/or  ring  current, 
precipitation  of  auroral  particles  in  an  increasing  longitudinal 
sector  of  the  midnight  auroral  oval,  and  finally  a poleward 
curvature  of  arcs  in  the  premidnight  area  that  resembles  the 
drift  pattern  in  the  polar  cap.  Minor  signatures  of  instability  or 
activity  are  followed  by  the  sharp  breakup. 

It  is  apparent  from  the  observations  presented  here  as  well 
as  from  observations  of  Subbarao  and  Rostoker  [1973],  Hones 
et  al.  [1973],  Akasofu  et  al.  [1973],  Rossberg  [1974],  and  Aka- 
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Fig  20  High-energy  eleciron  fluxes  averaged  over  pilch  angle  The  low-latitude  limit  of  the  position  of  the  boundary 
between  open  and  closed  field  lines  lies  between  invariant  latitudes  67.9°  and  68.4°,  indicated  by  hatching.  The  arrow  on  the 
absossa  labeled  01  IS  indicates  the  beginning  of  minute  IS  of  hour  01  (UT). 
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sutu  Iha  some  subsiorms  hase  .1  (ir.nliial  ami  easiix 

defined  growth  or  deselopmeni  phase,  while  others  ilo  not 
Ohsiousls.  ther  - is  no  elcar  one-to-one  relation  between  the 
growth  phenomi  non  ■•'■I  the  suhstorm  It  is  in  agreement  wnh 
current  substorn.  ineories  that  a /r,-negatise-saused  magneto- 


spherie  eonveelion.  if  it  is  eoiuinued  lor  a siillieientlx  long 
period,  mat  result  in  substorm  aetitiit.  as  in  the  ease  studied 
here  lliiweter.  the  etulenee  diseussed  abitte  seems  tii  argue 
against  the  use  vsf  the  term  growth  phase  when  it  is  applied 
with  the  same  weight  as  the  term  ‘expansive  phase'  in  describ- 
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ing  a definite  phi  ve  '■if  every  substorm.  The  term  may,  how- 
ever, be  suitable  in  describing  a magnetospheric  situation 
which  may  lead  to  the  triggering  of  a substorm. 

Ionospheric  Electric  Fields 

In  a study  of  the  ionospheric  electric  field  by  means  of 
balloon-borne  detectors,  Mozer  [1971]  found  a regular  behav- 
ior of  the  electric  field  in  the  auroral  zone.  Averaging  over  19 
substorms,  he  found  that  about  1 hour  before  the  onset  of  the 
expansive  phase  the  westward  component  of  the  electric  field 
increased  from  zero  and  reached  a maximum  value  at  the  time 
of  breakup.  The  field  remained  at  this  level  throughout  the 
expansion  phase  and  then  began  a return  toward  zero.  The 
southward  component  remained  at  a zero  level  until  it  sud- 
denly increased  at  the  onset  of  the  expansion.  Our  observation 
that  the  equatorward  drift  of  the  auroral  arc  appears  during 
the  growth  phase,  continues  through  the  expansive  phase,  and 
then  suddenly  stops  is  in  agreement  with  Mozer’s  observation 
of  the  duration  of  the  positive  excursion  of  the  westward 
component  of  the  electric  field.  Our  observation  of  the  limita- 
tion of  the  equatorward  drift  to  the  evening  side  of  the  sub- 
storm center  appears,  however,  to  be  new.  It  is  not  obvious 
from  Mozer's  balloon  measurements.  His  conclusions  are 
based  on  average  values  without  information  about  the  loca- 
tion of  the  balloons  in  relation  to  the  substorm. 

Particle  Precipitation  and  Relationship 
to  the  Magnetosphere 

The  particle  precipitation  associated  with  this  auroral  dis- 
play has  been  studied  both  by  direct  measurement  and  by 
observation  of  the  drift  of  cloudlike  auroral  forms.  Most  of  the 
electrons  responsible  for  the  luminosity  in  the  premidnight 
sector  before  breakup  possessed  energies  of  1-5  keV.  The 
cloudlike  form  observed  drifting  from  Tromsd  to  Murmansk 
was  produced  by  electrons  of  ai40-keV  energy.  This  energy 
determination,  based  on  the  assumption  that  the  electrons 
were  gradient  drifting,  is  in  agreement  with  the  photographic 
observations  of  Stormer  [\9S5\,  which  give  a mean  height  of  94 
km  for  ordinary  cloudlike  auroras  over  Norway.  The  ioniza- 
tion production  maximum  for  40-keV  electrons  has  been  cal- 
culated by  Rees  [1963]  to  be  92  km.  Our  observations  of 
diffuse  aurora  indicate  that  such  electrons  appeared  nearly 
simultaneously  all  along  the  oval  from  midnight  to  shortly 
before  0400  geomagnetic  time  and  that  the  drift  of  the  cloud- 
like aurora  east  of  TromsrI  is  consistent  with  gradient  drift  of 
40-keV  electrons  at  auroral  latitudes. 

This  observation  agrees  well  with  the  model  for  precipi- 
tation into  the  auroral  zone  of  electrons  in  the  >30-keV  range 
presented  by  Sletten  et  al.  [1971]  on  the  basis  of  observations 
by  several  authors  (see  their  list  of  references).  In  their  model, 
electrons  within  the  loss  cone  are  precipitated  in  an  elongated 
region  following  the  auroral  oval  from  geomagnetic  midnight 
toward  dawn.  Electrons  with  other  pitch  angles  drift  eastward, 
following  the  L shells  from  this  source  region,  and  are  gradu- 
ally precipitated  by  pitch  angle  scattering  or  other  mecha- 
nisms. This  drifting  electron  population  may  be  responsible 
for  the  mantle  aurora  [Hoffman  and  Burch,  1973]. 

Thus  our  observations  support  a mechanism  which  would 
inject  40-keV  electrons  into  the  auroral  oval  during  the  ex- 
pansive phase  of  the  substorm.  The  observation  by  the  Vela 
satellites  [Hones  et  a!.,  1973]  of  an  enhancement  in  the  E > 45- 
keV  electron  flux  at  the  time  of  rapid  thinning  of  the  plasma 
sheet  at  the  satellite  suggests  that  these  electrons  may  be  of  the 
same  population  as  that  responsible  for  the  production  of  the 


diffuse  cloudlike  auroras  which  appear  at  auroral  latitudes 
during  breakup. 

Our  present  Isis  I electron  observations  made  before  auroral 
breakup  agree  well  with  the  growth  phase  morphological  pat- 
tern of  Hoffman  and  Burch  [1973].  The  term  growth  phase  was 
defined  by  them  as  the  time  between  the  southward  turning  of 
the  IMF  and  the  appearance  of  auroral  breakup.  The  studies 
of  Akasofu  et  al.  [1973]  and  Akasofu  [1974a]  would  have 
categorized  the  present  observation  as  the  ‘plasma  sheet  quiet 
arc’  of  a ‘growth  phase  pattern.’  Winningham  et  al.  [1975], 
using  Isis  I and  2 particle  data,  have  amplified  the  patterns  of 
Hoffman  and  Burch  [1973]  but  do  not  include  a pattern  repre- 
senting a preexpansive  phase,  IMF  south  condition  (since  this 
condition  does  not  always  initiate  a complete  substorm  se- 
quence). They  do  state,  however,  that  an  observation  of  the 
kind  we  report  would  be  categorized  as  representing  their 
‘boundary  plasma  sheet'  (BPS)  region  during  a period  for 
which  the  IMF  is  directed  southward  (see  their  Figure  19). 
This  BPS  region  is  always  located  poleward  of  the  region  of 
diffuse  aurora  and  often  contains  discrete  arcs. 

As  we  have  illustrated  in  Figure  19,  the  highest  fluxes  of 
protons  of  about  4-keV  energy  were  measured  within  a region 
equatorward  of  the  electron  arc.  A similar  latitudinal  profile, 
characteristic  of  the  premidnight  hours,  has  been  obtained 
statistically  with  airborne  [Eather  and  Mende,  1971]  and 
ground-based  [Fukunishi,  1975]  Hd  observations.  The  study  of 
Fukunishi  [1975]  clearly  places  the  belt  of  proton  precipitation 
equatorward  of  discrete  auroral  arcs  during  times  which  he 
dt,.  ICS  as  the  ‘quiet’  phase  and  the  ‘prebreakup’  phase.  Our 
particle  observations  show  that  the  separation  can  occur  on  an 
instantaneous  basis,  a fact  which  must  be  included  in  any 
theory  which  seeks  to  explain  the  global  development  of  an 
auroral  substorm.  Isis  I crossed  the  arc  at  a local  time  of  20.5 
hours  and  a magnetic  local  time  of  20.8  hours.  We  also  note 
that  apparent  separation  between  the  region  of  maximum 
proton  precipitation  and  the  region  of  maximum  electron  pre- 
cipitation (the  arc)  does  not  completely  exclude  protons  from 
the  electron  precipitation  region;  it  only  excludes  protons 
within  the  energy  range  10  eV  to  1 1.6  keV  and  above  an  energy 
flux  level  of  a<IO'*  erg/cm*  s sr. 

The  observed  protons  were  most  likely  associated  with  the 
quiet  time  ring  current,  which,  according  to  the  San  Juan 
magnetogram  of  Figure  2,  intensified  as  the  substorm  sequence 
continued.  The  greatly  reduced  pioton  flux  at  the  position  of 
the  arc  before  breakup  implies  either  ( ) an  upward  field- 
aligned  electric  field  at  the  auroral  arc  or  (2)  that  the  mecha- 
nism which  energizes  electrons  and  protens  does  not  precipi- 
tate the  protons  efficiently.  Only  when  the  proton  population 
of  the  near-earth  plasma  sheet  reaches  a critical  intensity  are 
protons  precipitated  by  resonant  wave-particle  interactions  as 
suggested  by  Kennel  and  Petschek  [1966],  Adiabatic  compres- 
sion operates  on  protons  on  closed  field  lines  in  the  tail  region, 
but  the  dominant  mechanism  is  betatron  or  transverse  energi- 
zation which  operates  on  the  transverse  component  of  the 
proton  motion.  Our  observation,  using  the  high-energy  elec- 
tron data,  that  the  boundary  between  open  and  closed  field 
lines  lies  at  a latitude  at  least  as  high  as  67.9”  insures  that  the 
protons  we  observed  were  on  closed  magnetic  lines. 

Summary  and  Conclusions 

A comprehensive  study  of  the  major  substorm  of  February 
25,  1969,  at  0130  UT  and  several  weak  substorms  before  and 
after  this  event  has  resulted  in  the  following  observations. 

I.  The  preexpansive  phase  auroral  arc  extended  approxi- 
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mately  along  the  auroral  oval  from  17.5  to  05.5  MLT  and  was 
responsible  for  an  energy  input  rate  of  =>  3 X 10“  ergs/s  before 
the  breakup. 

2.  In  the  evening  sector,  from  the  earliest  observation  of 
auroral  light,  an  equatorward  drift  of  6 km/min  was  observed; 
this  drift  was  present  only  in  the  evening  sector  and  ceased 
after  the  onset  of  the  expansive  phase  at  the  time  of  maximum 
poleward  displacement  of  the  breakup  arcs. 

3.  In  the  evening-midnight  sector  the  increase  in  the  in- 
tensity of  the  preexpansive  phase  arc  accompanied  the  increase 
in  the  deviation  of  the  horizontal  magnetogram  component. 

4.  The  equivalent  currents  obtained  from  the  observed 
magnetic  deflections  were  of  the  DPZ  type  (twin  vortex  mode) 
before  the  onset  of  the  expansive  phase.  After  the  onset  the 
OPZ  currents  were  more  intense  by  a factor  of  2,  but  the 
dominating  current  during  this  period  was  an  intense  west- 
ward auroral  electrojet. 

5.  During  the  expansive  phase,  symmetrically  traveling 
disturbances  were  observed  propagating  eastward  and  west- 
ward away  from  the  midnight  sector.  The  propagation  stopped 
for  1-2  min  at  the  time  of  maximum  poleward  expansion  and 
then  continued  (at  least  in  the  morning  side)  with  the  original 
velocity. 

6.  F.lectrons  in  the  I-  to  5-keV  range  were  measured  by 
satellite  over  the  preexpansion  phase  arc.  They  were  observed 
near  and  probably  within  the  low-latitude  limit  of  the  bound- 
ary between  open  and  closed  field  iines. 

7.  Protons  of  2!4-keV  energy  were  measured  equatorward 
of  the  electron  arc;  their  intensity  was  1.6  X 10  ’ that  of  the 
electrons  within  the  arc,  and  they  were  definitely  on  closed 
field  lines. 

K.  An  injection  of  =^40-keV  electrons  as  determined  from 
the  observations  of  diffuse  cloudlike  aurora  occurred  during 
the  expansive  phase.  At  times  later  than  04  MLT  the  cloudlike 
aurora  associated  with  the  40-keV  electron  population  was 
ob.served  to  drift  eastward  toward  later  morning  hours. 

The  major  substorm  at  0130  UT  and  a smaller  substorm  at 
0600  UT  show  growth  phase  features  in  their  visual  and  mag- 
netic signatures.  However,  other  investigations  show  that 
many  substorms  have  occurred  for  which  there  is  no  classic 
growth  pattern  and  also  that  some  growth  phase  patterns  do 
not  always  result  in  a complete  substorm  sequence  [Kokubun  et 
al..  1977].  Accordingly,  we  emphasize  the  usefulness  of  the 
growth  phase  concept  in  describing  a magnetospheric  condi- 
tion which  is  initiated  by  a southward  turning  of  the  IMF  but 
stress  that  in  general  it  should  not  be  expected  to  bear  a one- 
to-one  correspondence  with  the  substorm  expansive  phase. 
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